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W. T. Hardgrave 


ABSTRACT 


This report gives a line by line description of sessions using 
the Set-Theoretic Data System (STDS) to interact with technical data 
bases. The data bases contain data from actual applications at NASA 
Langley Research Center. The report is meant to be a tutorial docu- 
ment that accompanies ,r Set Processing in a Network Environment" [l]. 


This paper was a result of work performed under NASA Grants NGR 47-102 001 
and NSG 1068 while the author was in residence at ICASE, NASA Langley 
Research Center, Hampton, Virginia 23665. 



1.0 Introductory Remarks 


This report is directed to the user who has a collection of data 
that he wants to query and/or update hut who has little, if any, 
knowledge of either set theory or information systems. In only a few eections do 
we make the tacit assumption that the reader has any experience in computing 
and these cases involve only a basic knowledge of FORTRAN, We hope to 
demonstrate, through the use of three scenarios, the power of a set- 
theoretic system in aiding the user to isolate and examine sets of data 
that may he of interest to him. The particular system used in this report 
is the Set-Theoretic Data System (STDS)[2] developed hy David Childs 
and supported hy Set Theoretic Information Systems (STIS) Incorporated, 

STDS was run (remotely from a teletype in Hampton, Virginia) on the 
IBM 360/67 at Wayne State University, Detroit, Michigan, using Michigan 
Terminal System (MTS) [ 3 ] developed at the University of Michigan, Since 
STDS is a research system rather than a production system, it is 
reasonable to expect it to have a number of deficiencies. Few of these 
can be classified as serious problems j most reside in the area of human 
factors and although revisions here are desirable, they are not absolutely 
necessary to solve the problems described herein. 

A standard strategy that pervades all of the scenarios, regardless 
of the data base is as follows. First, isolate the basic sets that can 
be requested by a single condition (e.g. velocity < 100 MPH). Then, 
using the set operations, isolate sets that have characteristics that are 
of particular interest to the interrogator. Whenever a new set is formed, 
the system declares the cardinality (i.e, the number of elements) of the 
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set. If the user desires more statistical information, he may request 
a summary of pertinent items that includes the maximums , minimums, means 
and standard deviations. If the user requires more specialized information 
about a given set, he may write a short FORTRAN routine that traverses 
a set, performs calculations and outputs the necessary values* An 
example of this usage is given in the VGH scenario. 
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2.0 An Introduction to STDS 


In this section, a brief introduction to STDS [2] is given. This 
consists of two subsections. Subsection 2.1 provides a description of 
the commands necessary and responses encountered when using STDS. In 
addition, other aspects of STDS that are common to all of the scenarios 
(e.g. timings) are dealt with in this subsection. 

Subsection 2.2 contains brief descriptions of some elementary STDS 
operations. Since STDS is primarily a research system maintained by 
STIS Corporation, it has a number of operations that are not particularly 
germane to the discussions considered herein. As a result, we list here 
only the operations that are relevant for a relational/set-theoretic 
system that might appear in a scientific user environment. These operations 
are listed in three categories: utility operations, restriction operations, 

and set operations. Subsection 2.2 is primarily meant to be used as a 
reference to accompany the examination of the scenarios. 

Finally, we offer a few opinions concerning the possible installation 
of STDS on CDC 6000 machines under KRONOS. STDS is written primarily in 
FORTRAN with some critical routines written in assembly code. The original 
STDS was implemented on an IBM 360/67 under MTS and a later implementation 
is available for IBM 360/370 machines under OS. Our estimate is that STDS 
could be installed under KRONOS in about three months. No special peripheral 
processor routines should be necessary. The areas that might cause diffi- 
culty are program arrangement and Input/Output (I/O) management. The IBM 
360/67 is a virtual memory machine and the original STDS implementation used 
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that feature to minimize manual segmentation of the program and data areas. 
Therefore some effort may be necessary in order to reduce the field length 
of STDS under KRONOS to an acceptable level. This will involve the use of 
buffer management and overlaying techniques. The I/O philosophies of IBM 
and CDC are radically different. As a result the conversion of STDS I/O 
to KRONOS may be tricky. However calls to CIO using the direct access 
capability should be sufficient. Higher level CDC supplied l/o routines 
should be avoided. 
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2.1 STDS Operations 


The STDS operations are listed below in three groups; utility operations 
restriction operations, and set operations. It is necessary to provide 
some informal definitions of the terms, elements, domains, and entries. 
Consider the relation 

A = {<1,2, 3, U>, <5,6,7 ,8>, <9,6,8,7>) 
which may also lie written 

A: 1 2 3 It 

5 6 7 8 

9 6 8 7 

Informally, the columns are equivalent to domains, the rows to entries. 

A relation is a set of entries. Each entry is an n-tuple containing n 
elements. Thus the relation may be represented as a table of elements 
with the columns corresponding to domains and the rows to entries. 

2.1.1 Utility Operations 
Name : GET 

Calling Sequences: 

STDS*: GET (C, 'DSN') 

FORTRAN: CALL GET (C,'DSN') 

Arguments: 

C is a set name 
’DSN’ is a data set name 
Description: 

A set stored off-line on disk, tape, datacell, etc. is brought 
into core as set C for active inquiry. The contents of 'DSN* 
must in fact be a set, not an array. Any 'DSN' for which GET 
is used must previously have been stored by PUT. 



Name : PUT 


Calling Sequences : 

STDS* : PUT (A, ’DSN* ) 

FORTRAN : CALL PUT ( A , ' DSN ' ) 

Arguments: 

A is a set name 

r DSN' is a dataset name 


Description: 

Any set A that has been created during a run using set operations 
may be saved off-line in direct access storage by use of PUT. The 
command PUT’ must have been used to place a set in storage before 
that set can be retrieved by means of the command GET. 


Name : ENTER 

Calling Sequences : 

STDS*: ENTER ( A, DA, ND A, 'FMT' ) 

Arguments: 

A is the resulting set 
DA is the domain declaration for A 
NDA is the number of domains in A 
FMT is the I/O FORMAT for A 
Description: 

ENTER allows a set. A, to be constructed by specifying the 
individual entries. ENTER is terminated by an END OF FILE. 
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Name: LIST 


Calling Sequences: 

STDS*: LIST(A,I,J t , FMT f ) 

FORTRAN : CALL LIST (A, I , J , ' FMT * ) 

Arguments: 

A is the set to he listed 

I,J are respectively the I-th and J-th entries of A 
'FMT 1 is a string name specifying the output format 
Description: 

This operation allows listing subsets of a set A under a format 
'FMT', The subset is the I-th through the J-th entries of the 
set A. 

Name: LISTU 

Calling Sequences: 

STDS* : LISTU 

Arguments : 

None 

Description: 

This operation allows a summary of the characteristics of all sets 
in the universe to be printed. 
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Name: FREE' 


Calling Sequences: 

STDS*: FREE (A) 

FORTRAN : CALL FREE ( A ) 

Arguments : 

A is the set to he freed 
Description: 

FREE{a) deletes A from the universe and frees the memory of 
its contents. 

Name: SETFMT 

Calling Sequences : 

STDS* : SETFMT (DD, ’ FMT 1 ) 

FORTRAN : I = SETFMT ( DD t • FMT * ) 

Arguments : 

DD is the domain declaration 8 5 l6, 32 
’FMT* is a FORTRAN format specif i cat ion 

I is the index associated with the specified format 

Description: 

SETFMT associates an index , I s with a format so that the format 
may he specified via '/I/'. 


OKIGUVAL ?App ra 

0i - 
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Name: SYSTEM 


Calling Sequences: 

STDS*: SYSTEM 

FORTRAN: CALL SYSTEM • 

Arguments: 

None 

Description: 

SYSTEM returns control to the command system in the host 
operating system. 
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Name: SUM 


Calling Sequences: 

STDS*: SUM ( A , INDEX , STAT , F ) 

FORTRAN: CALL SUM (A, INDEX, STAT, F) 

Arguments : 

A is the set to be summarized 
INDEX is the index set as described below 
STAT is the statistical set (NUMD=6) to be generated 
F is the floating point indicator; 0 implies integer 
arithmetic; 1 implies floating point arithmetic 

Description: 

SUM produces a statistical summary 0 f specified elements. 

Index Set : 

The index set should be DD = 32 NUMD = 3. Each triple <x,y,z^ 
in the index set should have the following form: 

X is an integer key. This key is necessary when a number of 
statistics are being gathered with a single index set, 

Y is the starting byte of the domain. For 32 bit values, the 
starting byte for a domain I can be calculated as follows: 

Y = 4*(l-l) + 1 

Z is the length in bytes of the domain. For 32 bit values, 

z = 4. 

Example: 

CALL SUM (A, INDEX, ST AT, 1 ) using 32-bit values and with 
INDEX = [<1,5>4> ,<2,9,4'>j would produce the sum, mean, 
minimum, maximum and standard deviation for domains 2 and 3 
(with keys 1 and 2 so the user can tell which is which). 
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Name : V$ 

Calling Sequence: 

STDS: V$($X,DD,V) 

Arguments: 

$X is the data mode. I implies integer, Z implies 

hexadecimal, A implies EDCDIC coded, and E implies 
floating point, 

DD is the domain declaration (usually 32) 

V is the name to he assigned to the constant 
Description: 

V$ provides the user with the capability to predefine 
constants to be used at a later time. After entering the 
V command, the system responds with a Then the user 

enters the constant. 

Example : 

V$($I,32,K0N) 

• 51 

The user enters the V$ command. The system responds with 
a The user enters the integer, ’57'. Thus the user 

creates the integer constant '57* with the alias KON. When- 
ever ICON is encountered in a command, STDS substitutes 57* 


12 



2.1.2 Restriction Operations 


Name: RSEG , RSLTA , RSLEA , RSGEA , RS GTA , RSNE 

Calling Sequences: 

STDS® : RS (l,A,V,C) 

FORTRAN: CALL RS (l,A,V,C) 

Arguments : 

I is a domain name 
A is the set to "be restricted 
V is a full word program variable 
C is the result set, a subset of A 
Description: 

The RS commands restrict A (forming C) as dictated by 

the domain I and the constant V. C is the subset of A 
containing entries of A whose Ith domain is in the 
specified relation (e.g. EQ denotes ’equal’) to V. 

Examples : 

For A = £<-10,20, 3> »< 0,5,10> ,< 6,7,1000>'j 

a) RSLTA ( 1 , A ,-2 ,C ) produces a set C = 1 < -10,20,3 

b) RSLEA(2,A,10,C) produces a set C = (<0,5,10 >,< 6,7 ,1000 >} 

c) RSGTA(2,A,10,C) produces a set C = (<-10,20,3 >} 

d) RSGEA(3,A,1000,C) produces a set C = }<6,7,1000>} 

e) RSNE(l,A,10,C) produces a set C = A 


ORIGINAL PAGE IS 
0E POOR QUALITY 
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Name : BRSA 


Calling Sequences: 

STDS*: BRSA(l,A,J,K,C) 

FORTRAN: CALL BRSA(l,A,J,K,C ) 

Arguments : 

I is a domain name 
A is the set to be restricted 
J,K are four byte program variables 
C is the result set, a subset of A 
Description: 

The set C is a subset of the set A such that the values of the 
entries in the I-th domain of C are greater than or equal to 
J and less than or equal to K, The set A is unchanged. 
Examples : 

if A = {<1,2,3,U>,<1,3,5,2>,< 2,3,6,7>} 

a) then BRSA(2,A,1 ,2,C) produces the set 
C = {<1,2, 3, U>} 

b) then BRSA(3*A,1+ ,10 ,C) produces the set 
C ={<1,3, 5, 2 >,<2,3,6,7>} 

c) and BRSA(h,A,U,10,C) produces the set 
C = {<1,2,3,^ >, < 2,3,6,7>} 
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Name: ALEG 


Calling Sequences: 

STDS*: ALEG(A,I,J,LEN,CL,CE,CG) 

FORTRAN: CALL ALEG(A,I,J,LEN,CL,CE,CG) ■ 

Arguments : 

A is the set to be partitioned 
I,J are domains 

LEN is the length of the domain in bytes (usually 4) 
CL,CE,CG are the output sets 
Description: 

The set A is partitioned into three sets, CL, CE and CG according 
to domains I and J. For each entry, K of A, if domain I is 
less than domain J then the Kth entry is placed in CL; if domain 
I is greater than domain J, then the Kth entry is placed in CG; 
and if they are equal, the Kth entry is placed in CE. 

Example : 

If A ={< 1,2, 3, 4 >,< 1,3, 5, 5 >,< 2,3, 7,6 >,< 5,3,2,7>} 

then ALEG (A, 3 ,4 ,4,CL,CE,CG) 

produces CL ={<1,2, 3, 4 >,< 5,3,2,7>) 

CE = {< 1,3, 5, 5 >} 

CG = {< 2, 3, 7, 6 >} 


original PAGE is 
a '' rooit QUALITY 
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Name : DMR 


Calling Sequences: 

STDS*: DMR(l,J,A,C) 

FORTRAN : CALL DMR(l,J,A,C) 


Arguments: 

I is a domain name for the set A 
J is the number of consecutive domains 
A is the given set 

C is the result set, the I-th domain of A 


Description: 

The set C is formed from the set A by taking the I-th through 
J-th domains of each entry of the set A. The set A is 
unchanged. 

Examples: 

For A = {<a,b,c>,< d,e,f > , < g,h,i>} 

DMR (1,2, A, C) produces 
C = {< a,b >,<d,e >,<g,h>} 
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Name : GETELM 


Calling Sequences: 

FORTRAN: CALL GETELM(l ,A, J,L,ELM) 

Argument s : 

I is the entry number 
A is the set name 

J is the first byte of the tuple to be extracted 
L is the number of bytes to be extracted (usually k) 
ELM is an array 
Description: 

GETELM extracts L bytes beginning at the Jth byte from 
the Ith entry of the set A. This information is trans- 
ferred to the user's array ELM. For 32-bit values (h bytes) 
the byte number may be calculated from the domain number by 
the formula BN=4*DN — 3 where BN is the byte number and 
DN is the domain number. 

Examples: (using 32 bit values) 

If A ={<1 ,2,3,1+ >, <1,3, 5,2 >,< 2,3,6,T >} 

a) Then GETELM (2»A S 5»^»ARR) produces 
ARR(l) = 3 

AKR(2) = 5 

b) Then GETELM (3,A,1,3,ARR) produces 
ARR(l) = 2 

ARR(2) = 3 
ARR(3) = 6 
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Name : SUBSET 


Calling Sequence: 

STDS: SUBSET(A,F,I,L,C) 

FORTRAN: CALL SUBSET(A,F,I,L,C ) 

Arguments: 

A is the input set 
C is the output set 
F is the first entry 
L is the last entry 
I is the increment 
Description: 

The Fth through the Lth elements in increments of I are 
copied from A thus forming C. 

Examples: 

For A = {<a,b,c >,<d,e 9 f >, < g,h,i>} 

SUBSET(A,1,2,3,C) produces 
C = {<a,b,c >, < g,h,i>} 
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2.1.3 Set Operations 


Name : UN 

Calling Sequences: 

STDS*: UN(AB,C.) 

FORTRAN: CALL UN(A,B,C) 

Arguments : 

A,B are names of given sets to be unioned together 
C is the name of the result set 
Description: 

Given two sets A and B, a new set C is produced containing each 
entry that is in either A or B. Note that if the same 

entry is in both A and B, only one representation of it 
appears in C. A and B are unchanged. 

Examples : 

If A = {< a,b,c >, < d,e,f > , < g,h,i >, < j ,k s l>} 
and 

B={< a,b f c >,<m,n,o >, < p,q,r >,< s,t,u>} 
then UN(A,B,C) produces 

C ={<a s b,c >,<d,e,f >,<g,h s i >,<j,k,l >,<m,n,o >,<p,q,r >,<s 
Note that < a,b,c > is in both A and B, and occurs but once in 
set C. 
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Name: IN 


Calling Sequences: 

STDS*: IN(A,B,C) 

FORTRAN: CALL IN(A,B,C) 

Arguments: 

A,B are the given sets to be intersected 
C is the result set 
Description: 

The set C contains all those entries that are entries of 
both the sets A and B. The sets A and B are unchanged. 
Examples : 

1. For A = {< a,b 4 c >, < d,e,f >, < g,h,i>} 
and 

B = {< a,b ,c > < g ,h,i > < j,k,l>} 

IN(A,B,C ) produces 

C = {< a,b,c> < g,h,i>} 


20 



Name : RL 


Calling Sequences: 

STDS* : RL(A,B*C ) 

FORTRAN : CALL RL(A,B,C) 

Arguments : 

A 3 B are the given sets 
C is the output set 
Description: 

The set C contains all those entries of the set A that are 
not entries of the set B. The sets A and B are unchanged* 
Examples : 

For A = {< a,h,c >, < d,e 4 f >} 
and 

B « {< a,b,c >, < g»h,i >, < J ,k,l>> 

a) RL(a,B.,C ) produces 

C = {< d,e,f>} because < d,e 4 f> is not in B 

b) RL(B,A,C) produces 

. C ={<g > h,i >,< j,k,l>} because <g,h,i > and > 

are not in A 


ORIGINAL PAGE IS 
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Name; SD 


Calling Sequences: 

STDS*: SD(A,B,C) 

FORTRAN: CALL SD(A,B,C) 

Arguments : 

A,B are the names of the given sets whose symmetric 
difference is to he found 
C is the name of the result set 
Description: 

Given two sets A and B, a set C is produced that is the set of 
entries present in either A or B hut not in Both. That is, 

C contains the entries that the two sets do not hold in 
common. A and B are unchanged. 

Examples : 

For A = {< a,b,c >, < d,e,f >, < g,h,i>} 
and 

B ={<a,h,c >,<d,e,f >,<x,y 9 z>} 

SD('A,B,C) produces 

C = {< g,h,i >, < x,y,z>} 
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Name; CARD 


Calling Sequences: 

FORTRAN: I = CARD (a) 

Arguments; 

A is the set 
Description: 

Returns the cardinality (i . e. the number of elements) of the 
set A in location I. 

Examples : 

If A ={<1,2,3,!+ >,<1,3, 5, 5 >,<2,3,7,6>} 
then I = CARD (A) produces 1=3 
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2.2 Entering; STDS 


The MTS command RUN RBLG:STDS* calls STDS into execution. MTS 
issues the EXECUTION BEGINS message just prior to transferring control 
to STDS. When STDS becomes active, it generates a heading and a request 
for the user to enter a universe name. The universe is the set of all 
sets known to STDS and is of little importance to the user. A list of 
the sets in the universe may be obtained by using the LISTU command. 

The user enters a four-character universe name of his selection; however, 
he must take care not to use that name for some other set during the session. 
Now the user is ready to proceed with his requests. 

At this point, the user would normally access some data base (i.e. set) 
that he has previously saved as a permanent file with the PUT command. 

He may do this with the GET command; for example, GETCAj'FE') finds the 
permanent file with the name FE and associates it with the local set 
name, A, Next, we discuss the two lines of output that follow the GET 
command, or any other set creation command. Consider some sample lines 
that might follow a GET command; the first, 

DD = 32 NUMD = 9 CARD = 420 

indicates that the domain declaration is 32, the number of domains in the 
set is 9, and the cardinality of the set is 420. The domain declaration 
of 32 implies that values in the set are stored as 32-bit entities. The 
number of domains equal to 9 implies that each element of the set is 
a 9-tuple; or more informally, the set has 9 columns of information. 

A cardinality of 420 implies that the set has 420 elements. Informally, 
that means that the set has 420 rows. However, it is important to 
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remember that theoretically the rows are unordered and no two rows would 
have identical values in all domains. 

The second line printed after a set creation, enclosed in brackets 
as shown below, contains timing information. All times are in seconds. 
[0.005, 0.05 / 0.136, 0.90] 

The two numbers preceding the slash indicate STDS operation times; 
first CPU time and then elapsed (i.e. real) time. The numbers following 
the slash are total times including MTS overhead; first CPU time and 
then elapsed time. They will, of course, vary with the load on the 
computing system. 

This same GET command used 0.005 CPU seconds and was completed in 
0.05 seconds of real time (i.e. wall clock time). However, the CPU 
time used including MTS overhead (for virtual memory paging, etc.) was 
0.136 seconds and the total real time from receipt of carriage return 
to the issuance of the first character of output was 0.90 seconds. 
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3 . 0 Finite Element Scenario 


The following scenario describes a situation that might arise in a 
computer aided aircraft design environment using a set theoretic informa- 


tion system (in this case STDS). Assume the designer is working with a 
finite element model of a fuselage, a portion of which is shown in Figure 3-1. 
Furthermore, assume that the model has been run through a suitable finite 
element program (in this case, SKA? ) applying suitable loading and boundary 
conditions as shown in Figure 3-2. These input conditions may be repre- 
sented as relations (because of their tabular nature) and maintained under 


STDS. We have not chosen to do that in this scenario because our original 
data was output from (rather than input to) the finite element program SNAP. 
This results in a data base of forces and moments that the designer wishes 
to examine in order to evaluate the design. The set-theoretic data base 
examined in this session has nine domains with 250 entries. A portion of 
the data base is shown below where FX, FY, FZ represent forces and MX, 

MY, MZ represent moments. 


FE 

TYPE 

MODE 

FX 

FY 

FZ 

MX 

MY 

MZ 

1 

Beam 

1 

78. 

-71. 

0. 

0. 

0. 

5980. 

1 

Beam 

2 

-78. 

71. 

0. 

0. 

0. 

-964o. 

1 

Quad 

1 

-78. 

ro 

CD 

VO 

2l+l. 

-0. 

-0. 

-0. 

1 

Quad 

2 

-160. 

599. 

-696. 

-0. 

-0. 

-0. 

1 

Quad 

9 

136. 

-509. 

825. 

-0. 

-0. 

-0. 


2 6 
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BOUNDARY CONDITIONS 


NODE TX TY TZ RX RY RZ 


10 1 

9 0 1 

15 1 1 

17 0 1 

25 0 1 

33 0 1 

39 1 1 

1+10 1 

1+9 0 1 

57 0 1 


0 - FREE TO MOVE 


0 10 1 
0 10 1 
1111 
0 10 1 
0 10 1 
0 10 1 
0 111 
0 10 1 
0 10 1 
0 10 1 


1 - FIXED 


LOAD CONDITIONS 


NODE 

DIRECTION 

MAGNITUDE 

57 

3 

10000. 

58 

3 

8666. 

59 

3 

5000. 

6o 

1 

-10000. 

6l 

3 

-5000. 

62 

3 

-8666. 

63 ■ 

3 

-10000. 


Figure 3-2 



After initializing STDS, as shown in Figure 3-3, we are ready to 
obtain access to the finite element data base (called FE). This is 
shown in Figure 3-k, The command GET(FE,'FE') retrieves the FE data 
base from the MTS permanent file system. The data base had been 
previously input to STDS and stored in set-theoretic format with the 
PUT command. In order to get a feeling for the data base, we would like 
to print a few entries on the terminal. This may be accomplished with the 
LIST command. However, it is necessary to provide a FORTRAN format to 
the list command j and since formats are typically long and c umb ersome to 
enter, we may declare it in a SETFMT command. STDS RESPONDS WITH FMT = 1 
indicating that for future reference we may access this format with /l/. 

We list the first few elements of FE using E-format, and then, having 
discovered the magnitudes involved, we define a new format using the F— 
specification and list again using the new format, /2/. 

In order to obtain statistical summaries for the overall behavior of 
forces and moments, using the SUM command, it is necessary to define an 
index set that specifies which domains are to be tabulated. We have pre- 
defined a suitable index set "FE- INDEX" and stored it in the MTS permanent 
file system. We may recall that set with GET(lNDEX, "FE-INDEX" ) , Follow- 
ing that we perform a summary on FE. The index set was set up to tabulate 
domains U, 5, 6, 7, 8, 9 (as listed under KEYS) representing FX, FY, FZ, MX, 

MY and MZ respectively. From this we notice immediately that the highs, 
lows, and standard deviations indicate that there may be doubtful areas 
in FX, FY, FZ, and MZ. 
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*'RUN RALG ? STDS ♦ 

# EXECUTION AEG I NS 

* + SET-THEORETIC INFORMATION MANAGEMENT INTERACTIVE INTERFACE 
C** STOS-I ** VERSJUM32- 1 7-741 


♦ ** + ** *+ * |<*t * + + + #; +■ +** *** ** + 

♦ENTER 4-CHARACTER UNIVERSE NAME + 
SUN IV 

♦ UNIVERSE UN I V HAS SEEN CREATED * 

**+ ***■ +*+ * + * jc** *+ + * + * *** ** + 


"TAU-ON ’ 

C CORE -MAX ( 3 -PAGE ) : 16 16 161 


Figure 3-3 


* + 
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?GET( FE, 'FE' ) 


00=32 NUM0= 9 CARO= 

420 







C 0.005. 0.05/ 0 

.136. 0.903 







? SSTFMT C 32. • ( I 4. IX. A 4 . [ 4. 6 C 1 X. ES . 1 > > * 

) 






FMT= 1 








?LiST(FE J l^ 5 > >/!/* ) 







« 

1 BEAM 1 0.8E 02 

-0.7E 02 0.0 


0*0 


0.0 


0 • 6 E 04 

1 BEAM 2 -3.3E 02 

0 . 7 E 32 0.0 


0*0 


0.0 


-0. IE 05 

1 QUAD 1 -0.3E 02 

0.3E 03 3.2E 

03 

• 0*0 


- 0.0 


-0.(3 

1 QUAD 2 -0.2E 03 

0.6E 03 -0.7E 

03 

- 0.0 


-3.0 


- 0.0 

1 QUAD 9 0.1E 03 

-0.5E 03 0.SE 

03 

- 0.0 


- 0.0 


- 0.0 

? SETFMT ( 32. ’ ( 14. IX. A4. 

I 4 > 6 t 1X.F3.0) > ' 

> 






FMT= 2 








?LISTCFE,l,5,'/2/’) 








1 BEAM 1 73 . 

-7 I . 

0 . 


0 . 


0 . 

5930. 

1 BEAM 2 -73. 

7 1 . 

0 * 


0 . 


0 . 

-9643. 

t QUAD 1 -73* 

269. 241. 


- 0 . 


- 0 * 

- 0 . 

1 QUAD 2 - 1 60 « 

599. -696. 


- 0 . 


- 0 . 

- 0 . 

1 QUAD 9 136- 

-509. 325. 


- 0 . 

ft 


- 0 . 

- 0 . 


?G£T(INDEXi 'FE-INDEX' ) 
l)D = 32 NUMO= 3 CA« 0 = 6 

t 0.035. 0.04/ 0.145. 0.603 


?SUM( FE, INDEX, SI, l) 


**+ KEf 

SUM 


MEAN 

LOW 

HIGH 

STO 

4 

- . 1 42333E 

32 

-0.03 

-12600.00 

16433.00 

2932.33 

S 

0.879335E 

02 

0.21 

-47000.00 

61200.00 

7755.78 

6 

0 • 562503E 

30 

0.00 

-91800.00 

92933.00 

13156.59 

7 

3.3 


3.0 

- 0.00 

0.00 

0.00 

3 

0.0 


0.0 

- 0.00 

3.00 

0.00 

9 --374224E 06 

l):>=32 Ml 1M0= 6 CARO= 

-39 1.01 
6 

-234033.33 

323333-03 

42365.67 


C 0-525. 5.34/ 0.756. 44.953 


Figure 5~^ 


51 



Assume for this scenario that moments greater than 15000 and forces 
greater than 1500 are considered excessive and the designer wished to 
isolate cases where those conditions occurred. As shown in Figure 3-5 
we decide to isolate the positive moments in Z. We do this by using the 
RSGTA command (Restrict-Greater-than-Arithmetic ) and forming a new set, 

POSRZ. A statistical summary reveals that the mean positive MZ is 49028.05 
considerably above the tolerable level of 15000.0. In order to determine 
exactly which entries are excessive, we isolate those that are equal to 
or exceed the threshold with the BRSA command (Between Restrict - Arithmetic). 
This command forms a new set EXCESSMZ having those entries whose values in domain 
9 are between (and including) 15000.0 and 323000.0. STDS responds that 
there are 28 of those - more than we would like to examine in detail. 

Therefore, we isolate the worst ones in EX2. We find there are only 3 between 
200000.0 and 323000.0; so we list them and we find that they involve nodes 
T, 15 and 6l. 

In Figure 3-6, we decide to isolate the finite elements that connect to 
nodes 7» 15j and 6l; the worst offenders in positive moments in Z. Using the 
RSEQ command (Restrict Equal) on domain 3, we form the sets H0DE7, R0DE15 
and N0.DE61 and list their contents. 

We perform a statistical summary on N0DE61 in order to check that 
i.he Forces and moments about node 6l sum to zero and -5000.0 (see figure 3.2). 

We notice that the sums for FY and FZ are -3.0 and -4998.0 respectively. The 
dt signer may wish to examine his loads and boundary conditions to determine 
If LI icse discrepancies and justified. 
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7RSGTAC 9. FE, 0.0,POSMZ) 

0D=32 N'JMD= 9 CARO- 53 

C 0.025. 0.03/ 0.162. 1.163 


? SUM C POSMZ » IN0EX.S2. I > 


*** KE/ 

A 

5 

6 

7 

8 


SUM 

0.131670E 05 
0. 102513E 05 

I40380E-99 
0.32S261E- 16 
-. 140238E-33 


-MEAN 

342.77 

193.42 

-g.00 

0.00 

- 0.03 


0.259349E 37 49023*05 


00=32 NUMDs 6 CARO- 


C‘ 0.071. 


3.09/ 0.236* 


LON 

-9623.00 
- 16633.03 
-0.03 
0.33 
-0.33 
107.03 


43. 131 


HIGH 
9620*03 
1 6603*33 
0.00 
0.00 
3*3 

323300.00 


STO ** + 
2545.64 
4446.21 
3 . 30 
0.33 
0.00 
65361*26 


? 8RSA( 9. P0SMZ. 1 5033.3. 323000 • 0 . EXCESSMZ > 
00=32 NUMO= 9 CARO- 23 

[ 0.004. 0.00/ 0 . 092. 3*203 

?dRSA< 9 * EXCESSMZ *23 3 3 30* 3* 323033.0. EX 2 > 
00=32 NUMOs 9 CARO= 3 

C 0.033. 0*. 00/ 0.091. 0.973 


? L I ST ( EX2. 1 .3. ’ /2/* > 

6 SEAM 7 -4610. 
12 BEAM 15 -9620. 
46 BEAM 61 5210* 


13103. 

-16600. 

7080. 


0 . 

0 . 

0 . 


0. 0. 3233 3 3 * 

0. 0. 225300. 

0. 0. 2140M0. 


Figure 3-5 
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? RSEQ < 3 > 7 s ? 7>N0DE7 ) 

QD = 32 NUMO = 9 CARD* 3 

C 0*012.. 3.01/ 0.1 44. 0.40] 


?list<:node 7 . i * 3 . ' />/ » > 

6 SEAM 7 -4610. 

6 QUAD 7 4*6 1 0 * 

. 99 REAM 7 0. 


13100. 0. 

17230. 13000. 

0. -13000. 


?RSEQ(3> FE,15,NODE15) 

DD=32 NUMD= 9 CARD= 5 

C 0.012. 0.01/ 0 • 143. 1.72] 


?Li ST C NODE 1 5 » 1 , S. '/a/' > 


6 

QUAD 

15 

710300. 

-37433 

12 

BEAM 

1 5 

’ -9620. 

-16600 

12 

QUAD 

15 

-12600. 

-470(30 

99 

BEAM 

15 

0. 

0 

133 

beam 

15 

0. 

0 


32900. 

0 . 

-37600. 
13000. 
-8230 • 


? RSEQ < 3 ,FE,6l,NODE6l) 
i)D = 32 NUMO = 9 CARD = 5 , 

C 0*011. 0.01/ 0.144. 0*33] 


?L I ST C N0DE6 1 . 1 

40 QIJAD 6 1 

41 QUAD 61 

46 REAM 61 

47 BEAM 61 
91 REAM 61 


5. ’ / 2 / • > 

8200. 2200. 

-8793. -8790. 

5210. 7080. 

-4620. -493 . 

0 . 0 . 


392. 

-7620. 

0 . 

3- 

1730. 


? SUM<. NGDE6 1 » 1 NDEX* S» 1 

*rt KF / 


4 

5 

6 

7 

8 
9 

0D“32 

C 0 . 


3.0 

-•300000E 31 
- • 49 9803 E 34 
0.0 
0.0 
0.0 

6 CARDs 

018. 0.05/ 


0.3 
— 0.60 
-999.63 
0.0 
3.0 
0.3 
6 


0.224. 4 > 


-3790.00 
-8790.00 
-7620.00 
0 . 0 
0 , (3 

214300.00 
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0* 0. 323030. 

— 0 . - 0 . — 0 * 

0 * 0 . 0 . 


- 0 . 

0 . 

- 0 . 

0 . 

0 . 


— 0 * - 0 • 

0. 225030. 

-0 • -0 * 

0 . 0 . 

0 . 0 . 


- 0 . 
-0 • 
0 . 
0 . 
0 . 


-0 • -0 . 
-0 • -3 • 

0. 214003. 

0 • -21 4000 » 

0 . 0 . 


n I Gn 

3200.00 

7330.00 
1 73 3.30 

0 . O 

• 0 . 0 

2 1 4003.00 


S ii> r r 
6212.77 
5147.30 
3372.15 

- 0 • 0 

- 0.0 

135345.44 


Figure 5-6 
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In Figure 3-7 » we isolate the forces in X, Y» and Z that exceed the 
positive threshold point of 1500.0. Having isolated these in BIGFX, BIGFY and 
BIGFZ, we would like to find those entries where the threshold is exceeded in 
all three dimensions. This can he done by intersecting BIGFX and BIGFY 
forming BIGFXY and then intersecting BIGFZY with BIGFXY forming BIGFXYZ. We 
discover that only 12 entries exceed in all three dimensions and we list those 
entries. 

Finally in Figure 3-8 some information about timings and costs 
is provided. The total time spent by the designer in querying the data base 
was less than 20 minutes (as shown by elapsed time). The CPU time on a IBM 
360/67 was less than 20 seconds and the cost of the job was less than $4.00. 
Note that the dollar amounts given here are for a university environment and 
would be somewhat higher in a commercial environment. 
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7BRSAC4* FE, 1500 . 0, 16U00 . 0 , BIGFX ) 

D0=32 NUMQ= 9 CARD* 53 

C 0.017* 0.02/ 3.154* 1.111 

?3RSAC5* FE, 1500 . 0, 61200 . 0, BIGFY ) 

DD=32 NUMD= 9 CARD* 55 

C 0.017. 0.03/ 0.152* 0.223 

? URSAC 6» FE, 1500 . 0, 92900 . 0, BIGFZ ) 

DD = 32 NUMD = 9 CARD= 103 1 

£ 0.019* 0.02/ 0.154* 0*183 

? IN( 8 I GFX* 81 GFY * 3 1GFX / > 

DD = 32 NU:4Da 9 CARD= 33 

£ 0.014. 0.04/ 0.121* 0.303 


?IN(QIGFZ.BIGFXr*3IGFXrZ) 

D0=32 NUMD= 9 CARD= 12 

C 0.014* 0.02/ 3.116* 0.173 


7LISTC3IGFXVZ* 1*15* • /2/ * ) 


5 

QUAD 

13 

5890. 

5890- 

2493. 

-0* 

-0. 

-0. 

6 

QUAD 

7 

' 4610. 

17200. 

13003. 

-0. 

“0 » 

-0. 

1 1 

QUAD 

21 

6090. 

639 0. 

9283. 

-0. 

-0. 

-0 « 

12 

QUAD 

22 

9020 . 

33700. 

25200. 

-0. 

-0. 

-3. 

1 7 

QUAD 

29 

5800 « 

5330. 

21300. 

-3. 

-3. 

-0. 

18 

QUAD 

30 

7300. 

27300. 

20600 « 

-0. 

-0. 

-3. 

23 

QUAD 

37 

3370. 

3370. 

28300. 

-0. 

-0. 

- 0 * 

23 

QUAD 

33 

2040. 

2040. 

31100. 

-0. 

“0 ■ 

- 0 . 

30 

QUAD 

39 

16400. 

61200. 

92900. 

. -0. 

-0. 

*0 * 

36 

QUAD 

47 

3190. 

1 1930. 

46500. 

-0. 

-3. 

-0. 

41 

QUAD 

54 

7960 . 

7963. 

892 3 . 

-0. 

-0. 

-0. 

42 

QUAD 

55 

1 59 0 . 

5943- 

19300. 

-0. 

-3. 

-0. 


Fig-ure 3—7 
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7STSTEM 



fSSIGNOFF 



* OFF AT 01 J 38 . 35 

39-11-74 


* ELAPSED TIME 

13.657 MIN. 

S.91 

♦ CPU TIME USED 

19. 29 4 SEC. 

S I . 7 4 

'* CPU STOR VMI 

13.783 PAGE-MI N» 

SI .32 

♦WAIT STOR VMI 

33-974 PAGE-HR. 


♦APPROX. COST OF 

THIS RUN IS 53- 

67 

♦ DISK STORAGE 

13423.191 PAGE-HR. 

S3. 44 

♦ APPROX. REMAINING BALANCE I S86. 

04 


s 


Figure 3—8 



^.0 Dayfile Scenario 


In this scenario, we consider a data base consisting of dayfile data from 
the computer complex at NASA Langley Research Center. The original data, 
taken from the Dayfile Summary Tape for 01 July 197^, is organized by Jobs. 

Each entry in the data base represents a Job and each Job has over 50 variables. 

A data base, called DATA3, containing most of these variables is available; 

however, it is not used in this scenario because it is very cumbersome to 

use at an interactive terminal having only 70 or 80 columns. As a result 

a data base, called NEWDAY, with all of the Jobs contained in DATA3 but 

with fewer variables (domains) for each Job was created. The 

data base, NEWDAY, has domains as shown in figure h-1. Only THRU, CPTIME, 

FL, OS CALLS, LINES, ESTIME, TAPE and CRT will be referenced in this 
scenario. 

The first segment of the interactive session is shown in Figure 4-2. 

STDS is called with the RUN RBLG:STDS* command. We gain access to the pre— 
piously created set NEWDAY (which is our entire data base) by using the 
GET (NEWDAY, ’NEWDAY') command. The cardinality of NEWDAY is given as 
12^+3. We also use the GET command to gain access to two index sets, FPX 
and INX, which will be used in conjunction with the SUM command to generate 
statistical summaries. An SETFMT command is used to set up a format /l/ 
for use in the LIST command. Finally, the LIST command is used to print 
the first three elements of the data base as duplicated below with domains 
(column) headings. 
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NEWDAY DOMAINS 


DOMAIN 

NAME 

DESCRIPTION 

1-2 

JOBNO 

Jot Identifier (7 characters) 

3 

THRU 

Throughput - Total time (in minutes) that job is in 
the system (Integer) 

k 

CPTIME 

Total CPU time used in minutes (floating point) 

5 

CMTIME 

Residence time (in minutes) in central memory not in- 
cluding ROLLOUT time (floating point) 

6 

PPTIME 

PPU time used in minutes (floating point) 

7 

FL 

Field length from job card (octal) 

8 

OSCALLS 

Total number of 0/S calls (floating point) 

9 

COST 

Cost in dollars (integer) 

10 

LINES 

Number of lines of output (integer) 

11 

ESTIME 

Estimated time limit - job card (floating point) 

12 

CP 

Machine: A, B, C, D, or 1 

13 

TAPE 

Number of magnetic tape assignments (l character) 

lit 

CRT 

Number of CRT assignments (l character) 

15 

DC 

Number of data cell control cards in job (l character) 



Figure 1+-1 
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In Figure U-3 we call for statistical summaries of the entire 
data base with the SUM (NEWDAY, FPX, SI, l) command and the SUM (NEWDAY, INX, 
S0 a 0) command. In this version of STDS, it is necessary to ask for floating 
point summaries and integer summaries separately. Floating point domains 
are summarized with the FPX index set and a "1" in the last parameter posi- 
tion in the SUM command. Integer domains are summarized with INX index set 
and a "0" in the last parameter position of the SUM command. After the 
summaries, we calculate CP efficiency and average thruput for all Jobs. 

From other sources we determined that the total production time for 
all five machines was approximately 117 hours for 01 July 197 I*. We can 
calculate CP efficiency (as a percentage) by taking the sum of CPTIME 
(domain 4) converting it to hours, dividing by total production time (117 
hours) and multiplying by 100. The result is about 6l% CP efficiency. 

The CALC command is an MTS command (not an STDS command) that allows the 
teiminal user to evaluate expressions immediately. Average thruput for 
all jobs may be calculated by converting the average thruput (domain 3 ) 
to hours. The result is an average thruput of U.l hours. We use thruput 
as an approximation for turnaround since there are no figures available 
from vhich we can calculate turnaround. Turnaround is the time measured 
from the submission of a job to the time that the job is finished and 
is available for pickup by the user. Turnaround is equal to thruput 
plus the time necessary for manual handling of the Job. 



"HUN RALGjSTDS* 

* EXECUTION 3EGINS 




SET-THEORETIC INFORMATION MANAGEMENT INTERACTIVE INTERFACE 
C** STDS-i *♦ VERS:UMv32“ 1 7-743 




*** n*r +** *** * t ¥ * + * * + + *** **+' 

♦ ENTER 4 -CHARACTER UN I VERSE NAMF * 
SUN IV 

♦UNIVERSE UNIV HAS SEEN CREATED * 

* ♦ * * * * ('*♦ ,tt+ »»t + *• + # * + +»■* * + + 


‘'TAU-ON'' 

ICORE-.MAX ( 3-PAGE ) s 16 16 161 


7GETCNEWDAT, ' NEWOAf ' ) 

00=32 NUMD= 15 CARD= 1243 
t 0.004. 0.01/ 0.1R7» 0.29] 

?GET<FPX> ‘DA/-FP INDEX* > 

00=32 NUMOs 3 CARD= 5 

C 0.005. 0.03/ 0.133. 0.3U 

?GET< I NX » ‘DAT-INTINDEX* > 

00 = 32 N'JMD= 3 CARD= 4 

1 0.005. 0.01/ 0.133, 0.29] 1 


? SETFMT C 32 . * < l X »2A4. I 4 . F7 . 2 . F6 . 1 . F7 .2 • I 7. F5» ! » 1 4 . 1 6 » F5 *0. I X. 4A1 ) * ) 
FMT= I 


? LIST C NE/ZDAr' . 1 
CT62591 1 54 

EXP0A01 34 
EXP0A02 735 


.3. * / 1 / ' > 
41*40 140*4 
0*3! 3*4 

1 *06 4*2 


40*51 152000 

\ 0 • 3S 52000 

2*30 145000 


1 TOli . 

n * o 

f -T r»i-> 

1 r»n • 



* ^ A 

* • 

1 f 

66. 

5 

0 

1 . 

0301 

1 576 • 

19 

-6835 

5. 

i2332 


Figure 4-2 
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?SUMC NEWDAf » FPX » S 1 # 1 ) 
*** KElf sum 

A 0-425530E 34 

5 0 • l 9 1 B23E 05 

6 0* 725954E 04 

« 0.271834E 07 

I l 0* 176702E 05 

00=32 N 11*10= 6 CARO= 

i 1 * 535 » 2*19/ 1 . 


■ MEAN 

LOW 

3-42 

0*02 

15.43 

0*07 

'5-34 

0.07 

2136-92 

1 *00 

14*22 

0*10 


5 

» 35 *151 


HIGH 

142*71 

357*03 

351*64 

56985*03 

1296*33 


?SUM( NEWOAV. I NX* S3. 3 ) 


*** KET . SUM 

3 3 • 333233E 06 

7 0 • 903563E 03 

9 0-353630E 05 

10 3*41 0342E 37 

C3-32 NUM0= 6 CAR0= 
C 1 * 322* 2-15/ 


MEAN 

247-97 

72692.50 

23*45 

3293.81 

4 

1-538* 31-361 


LOW 
1 *00 

7000.00 

5*00 

140*03 


HIGH 

1425.30 

342033*03 

2336.03 

82551.00 


?SCOMMENT CALCULATE CP ERE I C I ENC i AND AVERAGE THflUPUT 

‘7'^CALC < t • 42553E+4 '40 > / 1 1 7 j * | 00 
*=60-6163391 163 


STD ♦** 
9*41 
33*01 
« o * j9 
4615*98 
67*01 


ST£> *** 
223*73 
37931.91 
126*67 
6692*17 


?5C4LC 247-97/63 
*=4-132333333333 


Pigure if— 3 
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In Figure k-k we have the interactive segment in which we isolate 
the jobs using one minute of CPTIME or less. We consider this to be a small 
demand for the CP resource. We isolate these jobs in a set called SMALLCP 
and notice that the cardinality of SMALLCP is We may then perform the 

statistical summaries on SMALLCP, We may calculate the percentage of jobs 
by using the cardinalities of SMALLCP and NEWDAY. The result is that about 
60% of the jobs are SMALLCP jobs. We then use the SUM of CPTIME (domain 4) 
for SMALLCP and NEWDAY to calculate percentage of CPTIME used. The result 
is that the SMALLCP jobs only use k.3% of the CP resources. Finally, we 
calculate the average thruput for SMALLCP Jobs by converting mean THRU 
(domain 3) to hours. The result is 3.2 hours. 

In the next segment of the session. Figure U— 5 we isolate Jobs using five 
minutes or less of CFTIME. We can classify those as small and medium length 
Jobs. We isolate these by using the RSLEA command in the set MEDCP. Rote 
that SMALLCP is a subset of MEDCP, The cardinality of MEDCP is 105 U, After 
asking for the floating point and integer summaries, we can calculate the 
job percentage, the CP percentage and the average thruput in the same manner 
as before. The results are that about 85 $ of the job load are MEDCP jobs. 

They account for only 2i% of the CP resources and their average turnaround 
is 3.7 hours. 

In the Figure h-6 segment, we isolate the jobs that do not require tape 
mounts. The command V$($A, 32, ZERO) sets up a constant that is type 
character (as opposed to integer) and has a value of 'O'. The domains 
TAPE, CRT, and DC were inadvertently input to STDS as type character 
rather than type integer; therefore it is necessary to retrieve it in 
character mode. The command RSEQ( 13, NEWDAY, ZERO, NOTAPES) isolates the 
no-tape jobs in NOTAPES. The cardinality is 806 . After the two standard 
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75C0MMENT ISOLATE JOdS WITH CPTIME 
?RSLEA< 4#NEW0Af# I .0#SMALLCP> 


00=32 NUMO* I 5 CARDs 74 3 

t I 40# 0.72/ 0.591# 4.71] 

? SUM ( SMALLCP# irPX # S 1 » 1 > 

KE/ SUM MEAN 

4 0.I83195E 03 0.24 

5 3 • 338022E 04 4.52 

6 0.I9I032E 04 • 2*55 

8 0 « 56662 1 E 06 757.51 

11 0.451563E 34 6*34 

D0=32 NUMU= 6 CARD* 5 


C 0.755# 1.92/ 0.971# 40.80] 


•LE. 1 .0 MINUTES 


LOW 

HIGH 

STU 

0.02 

0.93 

0.26 

0.07 

116.23 

3.04 

0.07 

54*15 

4,31 

1*03 

1 7447.0O 

1679-99 

0.10 

75.00 

l *3 . 49 


? SUM( SMALLCP# I NX. S3. 0) 


* + + KE/ sum 

3 0. 143376E 06 

7 0.431597E '33 

9 0.436133E 34 

10 0.133142E 07 

130=32 NUMi)= 6 CARO* 
C 0 . 643» 5*39/ 


MEAN 

192*35 

64384.63 

6.53 

1445.75 

4 


0.352. 33. 


LOW 

1 .00 

7000.00 

5.00 

140.03 


39] 


HI GM 
1425.33 
230000.00 
50 « 00 
23609.03 


STO * + * 
132.65 
32335.48 
4.22 
221 1 .43 


? '■’COMMENT CALCULATE %JOdS JCP AND AV. THRUPUT 


75CALC 100*743/1243 
■*' =60 . 17 6991 1 5044 


?$CALC 100*183.195/4255.3 
V =4*3051 01872953 


? •' C ALC 1 9 2 - 3 5 / 6 0 
* =3 *235333333333 


Figure L— L 
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75C0MMENT ISOLATE J03S WITH CPTIME .LE. 5.0 MINUTES 

? RSLEA ( A * NEVJOAT * 5 * 0 » MEOCP ) 

DD=32 NUMU- 15 CARO = 1054 

l 0 . 1 73 * 2.60' 0*641* 9.331 


?SU''U«£OCP*FPX»Sl* 1 > 


++* KEY 

SUM 

MEAN 

LOW 

H I GH 

STD 

A 


0.90139SE 03 

0.86 

3.02 

4.86 

1 .14 

5 


0 • 853943E 04 

3*15 

0.07 

357.03 

18.39 

6 


0.449352E 34 

4.26 

0.97 

351 .64 

12.35 

8 


0*1 4999 1 E 37 

1423.06 

l .00 

31043.03 

2752.34 

1 1 


0.1 3307 5E 05 

9.49 

0.10 

1296*33 

58.36 

D0=32 

NUMD = 6 CARD = 

5 




C 1 . 

074* 6.63/ 1 

•293* 47 

.42] 



? SUM ( MEDCP* INX*S3*0J 





*** KET 

SUM 

•MEAN 

LOW 

, HIGH 

STD 

3 


0*23421 IE 06 

222.21 

1.00 

1425.00 

206.30 

7 


0.729292E 08 

69192.75 

7 0 0 0 * 3 0 

240300.00 

35328-25 

9 


0*11 3230S 35 

10.74 

5.00 

965.00 

32.23 

10 


0.2S6065E 07 

2429.46 

1 40 .00 

27633.00 

3394.53 


DU=32 NUMO = 6 CARU = 4 

l 0*909* 2*58/ 1*107* 33.07] 


?f, COMMENT CALCULATE 2J06S* '.CP* AND AV. THRUPIJT 

76CALC 100+1054/1243 
#=34.79435116653 


#=21 .13287782295 


? -iCALC 222.21/60 
# =3*7335 


Figure 4—5 
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? ^COMMENT ISOLATE NOTATE JOdS 


7V5C 5A*32»ZERO> 

*0 


? RSE£ C I 3 , NE ‘tIO A If » ZERO . NO TAPES ) 
00=32 NUMD= IS CARO= 836 
C 0*057* 0.09 ' 0 • 49 I » 1 


401 


?SUMC NOT APES* EPX* S 1 * L > 


+ *r * KET 

SUM 

MEAN 

\ n m 
■! 

H I GrI 

QTn 

w i 

4 

0 . 227963E 04 

2*83 

0.02 

69.91 

7.49 

5 

0.101005E 05 

12*53 

0.07 

357.03 

29.30 

6 

0.319S57E 04 

3*97 

0.07 

351 .64 

13.35 

3 

0 . 855233E 06 

1061 .08 

1 .00 

35512.00 

3047.36 

1 1 

0.109764E 05 

13.62 

0* 10 

1296.30 

74.83 

DD=32 NUMOs 6 CAROs 

5 




C 0.811 

. 0.88/ 1 

.019* 36. 

781 



? SUM (NOT APES* INX»S0*0> 


\ 



+ + * KEY 

SUM 

MEAN 

LOW 

MI GM 

STD 

3 

0.156442E 06 

194.10 

1 .00 

1150.30 

172*34 

7 

0.577100E 08 

71603.33 

7000 . 00 

250000.00 

36971.63 

9 

0.169150E 05 

20.99 

5.90 

1 499 .00 

92.09 

13 

0.180654E 07 

2241 .36 

143* 00 

26965.30 

3376.33 


DD=32 NUMDs 6 CARO* 

C 0.697* 0.71/ 0.389* 


+ * * 


29.51 1 


? ^COMMENT CALCULATE SJ03S* %CP, AND A V. THRUPUT 

? SCALC 100*306/1243 
#=64.34312148023 


? AC ALC 100*2279.63/4255.3 
•7 = 53 *57 1 54607 t 9 5 

7SCALC 194.1/60 
7=3*235 


Figure b~6 



summaries are generated, we calculate the Job percentage, the CP percentage 
and the average thruput. The results are that about 65 % of the Jobs are 
no-tape Jobs. These require $h% of the CP time and have an average 
thruput of 3.2 hours. 

In the Figure 4-7 segment , we isolate the jobs requiring one or more 
tapes. This is easily done by "subtracting'' the set NOTAPES from the set 
NEWDAY and forming the set TAPES. The relative complement command RL 
performs this function. After the summaries are requested, we again calculate 
the Job and CP percentages and the average thruput. The approximate results 
are that 35% of the jobs are tape jobs. They require b6% of the CP resources 
and the average thruput is 5.9 hours. 

At this point we would like to isolate the jobs that fit into our 
general concept of a small job. The criteria given here is somewhat arbitrary, 
but it serves to illustrate the power of a set-theoretic information system 
for isolating classes of jobs that the interrogator wishes to examine. 

The criteria for a "small job" is as follows: 

CPTIME <_ 1.0 minute 

FL < 100000 words (octal) 

LINES < 2000 lines 

OS CALLS <_ 2000 0/S calls 

and no CRT's assigned. 

The Figure 4-8 segment demonstrates how this set called SMALLJ0B, is 
isolated and Figure 4-9 and Figure 4-10 shows the analysis of small no-tape 
jobs and small tape jobs respectively. In Figure 4-11 we isolate the 
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’^COMMENT ISOLATE TAPE J033 

?RL < NEMO AlT » NOTAPES # TAPES ) 

00=32 NUMD= 15 CARO* 437 
C 0.372# 1.57/ 0*535# 1.813 


?SUMCTAPES#EPX#Sl # t ) 


++*, KEY 

SUM 


MEAN 

LOW 

HIGH 

STJ 

4 

0. I97575E 

04 

4.52 

0.23 

142.71 

12.11 

5 

0.903263E 

04 

20.7 3 

0.13 

246.52 

• 30.57 

6 

0.406197E 

04 

9.30 

0.13 

94. S9 

I 1*73 

8 

3.18631 IE 

07 

4263.40 

2.00 

56935.00 

6063.27 

1 1 

0 * 669441 E 

04 

15.32 

0.20 

430 *00 

49.27 


UD=32 NUMO= 6 CARDs 5 


C 0.443# 0.49/ 0.643# 36*293 


? SUM (TAPES# I NX » S0 , 0 ) 


\ 



+ + + KEV 

SUM 

MEAN 

LOU 

HIGH 

STD 

3 

0.151783E 06 

347 .34 

3.03 

1425.00 

283*71 

7 

0.326474E 03 

74707.94 

7000.00 

342300.00 

39562.93 

9 

0.134450E 05 

42.21 

5.00 

2336.00 

172.35 

13 

0 • 229389E 07 

5249 . 1 7 

52 . 00 

32551 .00 

1 3024.36 

D0=32 NUMOs 6 CARQ = 

4 




C 0.331 * 0.40/ 0.572# 29*211 

75COMMENT CALCULATE .!J03S» %CP# AND Ay. THRUPUT' 

76CALC I 00*437 / 1 243 
* =35 . 1 5687351971 

? n if^a+ioTQ.it; / /,? r r . q 
*=46* 4302339365 

? -CALC 347 . 34 ''6.3 
*=5.739 


Figure k—'j 
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.? * COMMENT ISOLATE SMALL JOBS 

?RSLEA< 7 » NEWOAY, 1 00303 , SMALLFL > 

D0=32 iMU MD= 15 CARO= 993 

C 0 • 165. 3*41/ 0.610. 3.653 

? RSLEAC 1 0 , NEWDAY . 2000 » SMALL I iSiE ; 

00=32 NUMO = 15 CARl] = 770 

C 0.141» 0-35/ 0.578. 2.633 

? RSLEA ( 3 . NEW DAY , 2000 . 0 . SWALLOS > 

DD=32 NUMD= 15 CAR0= 951 

C 0.161. 0.21/ 0.602. 1.69] 

7RSEQC t 4.NEW0AY. ZERO. NOCRT) 

00=32 NUMO= 15 CARDs 1227 
C 0-077. 0.54/ 0.536. 2.123 

? I N ( SM ALLOA . SMALLFL . SM ALL JOB > 

00 = 32 NUMO= 15 CARl)= 658 
C 0.121. 1.33/ 0.262. 1.90] 

? IN( SMALL JOB. SMALLItNE. SMALL JOB > 

00 = 32 NUMOs |5 CAROs 522 
C 0.099. 0.94' 0.241. 1.151 

7 I N ( SMALL JOd , SMALLOS » SMALL JOB > 

Q0s32 NUMOs 15 CARO= 495 

T (A t Q . 4<J / n o/«t. r* /■-?-» 

- -i . j 4 * ^ t j 

7 1 N ( SMALL JOB , NO CRT . SMALL JOB ) 

00=32 NUMO= 15 CARO = 494 

C 4*106. 0.57/ 0.246. 0*763 


Figure 4-8 
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small field length Jobs in SMALLFL (card - 993), the small output jobs in 
SMALLINE (card a 770), the small 0/S call jobs in SMALLOS (card a 951), 
and the NO— CRT jobs in NOCRT (card — 1227). We may now form the inter- 
section 


SMALLJOB = SMALLCP (1 SMALLFL fl SMALLINE D SMALLOS fl NOCRT 

by performing four successive IN functions as shown. The (final) cardinality of 
SMALLJOB is U 9 U or about 39.7$ of the jobs. 

In the Figure k-9 segment, we isolate the small no-tape jobs in the 
set NOTAPESM. This is the class of jobs that one would expect to have the 
shortest turnaround. We perform the standard statistical summaries as shown. 

The HIGH figure for THRU (domain 3) indicates that turnaround can go over 
12 hours for even the shortest jobs. We calculate the job percentage, the CP 
percentage and the average thruput. The results are that the small no-tape 
jobs account for about 30# of the jobs and use only 1.3# of the CP resources. 

The average thruput is 2.1 hours. Note that the total CP time necessary to 
run these 376 jobs was less than one hour. One 6000 machine devoted to small 
jobs of this class during normal working hours would run between 1500 and 
3000 jobs in an eight hour period depending on the CP power and efficiency. 

In the Figure U-10 segment, we isolate the small, tape jobs by intersecting 
SMALLJOB with TAPES forming SMALTAPE (card - 118). The statistical summaries 
are requested. We calculate the Job percentage, the CP percentage, and the 
average thruput. The results are that almost 9 . 5 # of the jobs are 
tape jobs, they use almost 0.6# of the CP resources and the average thruput 
is k.k hours. In addition, the HIGH value of THRU (domain 3) indicates 
that turnaround for a small tape job can approach 2 l+ hours. 
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?^COMMENT isolate small jobs without tapes 


? INC small JOB. NOTAPES* NOTAPESM) 

00=32 NUMO= 15 CAHD= 37 A 
r 0.034, '3. 56' 0.21 2, 3.711 


7SUM C NOTAPESM, r PX , S 1 » l ) 


*** KEY SUM MEAN 

4 0.561680E 02 0.15 

5 0.100656E 04 2.63 

6 0.4987 17E 03 1 .33 

8 0.879460E 05 233.90 

11 0.1 93047 E 04 5-13 

00=32 NUMOs 6 C Af?0= 5 


t 0.382* 0.42/ 0.593, 36*303 


LOW HIGH STD *+* 

0.02 0.93 0.20 

0.07 109.55 6.56 

0.07 11*68 1.69 

1.00 1841.00 309.99 

0.10 62.90 9*68 


7SUMC NOT APESM* I NX* S3* 0 
*»* KEY SUM 

MEAN 

LOW 

HIGH 

3 

0.474293E 05 

126.14 

1 .00 

765.00 

7 

0.203439E 03 

53321 .53 

7000 * 00 

1 0000 0*00 

9 

0. 192133E 04 

5.1 1 

5.00 

9.00 

10 

0.133934E 06 

439.32 

140.00 

1973.00 

00=32 
£ 0- 

NUMD= 6 CARD= 
328* 0*33/ 0 

4 

*521* 29 

.221 



STD 
96.19 
14192. 30 
0.45 
571 .33 


7 '.COMMENT CALCULATE 2 JOB, 7.CP, AND AV. THKUPUT 

7 S CALC 100+376/1243 
V =30. 24939662 I 37 

? A CALC 100 + 56.168/4255.3 
*=l .319953939792 

7-CALC 126.14/63 
7=2*1 02333333333 


Figure 4-9 
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75COMMENT 


ISOLATE SMALL JOBS WITri TAPES 


? INC SM4LLJ03. TAPES. SMALTAPE > 

00=32 NUMDs 15 CARD = 118 

C 0.056. 0.64/ 0.176. 0.793 


7SUMCSMALTAPE.FPX.Sl . 1 > 


*+* KF9 SUM 

MEAN 

LOW 

HIGH 

STO 

4 0.251997E 02 

0.21 

0*23 

0.97 

0.23 

5 0.61 1853E 03 

5.19 

0.13 

39.50 

6.33 

6 0.272223E 03 

2.31 

0.13 

12.63 

2.20 

8 0 • 60937OE 05 

516*42 

2.30 

1913.00 

524.16 

11 0 » 339 1 95E 03 

7.11 

0-20 

50 * 00 

12*35 

00=32 NUMD= 6 CAR0 = 

5 




C 0.124. 0-15/ 3 • 

330. 36.033 



* 


7 SUM (SMALT APE. INX.: 

S3.0) 



r 


+ * + KEf 

SUM 


MEAN ' 

LOW 

HIGH 

STO + 

3 

0.3137 13E 

35 

265*8 6 

3.00 

1425.00 

229 « ! 3 

7 

3 • 6 4 3 6 48 E 

07 

54292*13 

7003 « 00 

I 00030.00 

18329.57 

9 

0 • 6470CQ0E 

03 

5-43 

5.00 

29.00 - 

2*52 

10 

0.67 1560E 

05 

569.12 

52.03 

1964.03 

548.31 


00=32 NUM0= 6 CARD* 4 

C 3.108. 0.11/ 0.299. 29.033 


7SC0MMENT CALCULATE /{JOBS. 7.CP, A V. THRUPUT 

? nCALC 100 + 118/1243 
#=9.49316 1 705551 

7 i-CALC 130+25-1997/4255.3 
# = » 592 1956143303 

75CALC 265.36/60 
*=4.431 


Figure 4-10 
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Figure 4-12 gives a manually created summary tabular form of the results 
obtained during the session. Figure 4-11 provides some information about the 
length and cost of the session. The total wall clock time for the session 
(as given by elapsed time) was less than an hour and the cost of the session 
(at university rates) was just over $10.00. 
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?SCOMM£tNT END OF SESSION 


7SVSTEM 
s' UN 

■^SSIGNOFF 


*OFF AT 33:51.53 
^ELAPSED time 
-CPU TIME USED 


39-1 3-74 
51*653 MItN. 
3S.775 SEC. 


^CPU STOP VMI 75-334 PAGE-MlN. 

'WAIT STOP VMI 115.319 PAGE-HR. 

DRUM RE ADS 1593 

APPROX. COST OF THIS RUN IS «!0-!! 

* D I SK STORAGE 3344.564 PAGE-HR. 

4 APPROX-. REMAINING BALANCE : •" SSI. 50 


02.53 

53.50 

54.33 


-il .37 


> 


Figure 4-11 



JOB CLASS 

# JOBS 

# CP 

AVG. THRUPUT 
(Hours) 

NEWDAY 

100.0# 

100.05? 

4.13 

SMALLCP 

60.185? 

4.35? 

3.21 

MEDCP 

84.79# 

21.185? 

3.70 

NOTAPES 

64.84# 

53.57# 

3.24 

TAPES 

35.16# 

46.435? 

5.79 

NOTAPESM 

30 . 255? 

1. 325? 

2.10 

SMALTAPE 

9.49# 

0.59# 

4.43 


CP Efficiency 

= 6l# 



Figure 4-12 
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5.0 VGH Scenario 


The first segment of the VGH interactive session is shown 
in Figure 5-l« We are in a position to retrieve the data base 
from mass storage and make it accessible to STDS. The GET 
command performs this function. Then, we list the first five 
elements in the data base. Ea,ch column in the listing represents 
one domain of the set. Domain number one is the aircraft type, 
domain two is the serial number, and so on as shown in Figure 5-2. 
Only the serial number, the gust indicator, the indicated air- 
speed, the altitude, and the acceleration a, re used in this 
scenario. 

For this scenario, we have prepared an appropriate index 
set in advance and stored it as a permanent file. We can re- 
trieve it with the command: GET( INDEX, 'VGH- INDEX' ) . The set, 
INDEX, provides for statistical summaries to be generated for 
the domains, IAS (8), Q(9), ALT(lO), P(ll), and ACC(l2). Now, 
we may generate a statistical summary of the entire data base 
(i. e. the set VGH) with the SUM command. 
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*run rrlgjstos* * 

♦EXECUTION 8E31NS ' • 

** SET-THEORETIC INFORMATION MANAGEMENT INTERACTIVE INTERFACE ** 
£ *•+ STDS- I ** VERS: UM32- 17-7 4] 


♦ * * ♦ * * * * * ♦ + * ***• tv* *** *** * * * 

♦ENTER 4-CHARACTER UNIVERSE NAME * 

: !JN IV 

♦UNIVERSE UN IV HAS HEEN CREATED ♦ 

♦ v+ * * * * * + * * + * + * * * * * + * * + + * * + 


'•TAU-ON'’ 

t CORE-MAX (8-PAGE) S 16 16 16] 

? SETFMT (32 , 1 <214* 1 5, 13# 12, 1 X»2 1 5> 5F3 .2 ) * ) 
FMT= I 

?GET(VGH, *HapT:VGH-SET' ) 


DD = 32 NUMD 

s 

12 CARDs 7263 






£ 0.005. 


0.01 / 

0.321 , 


0.75] 

. 




? LI ST (VGH, 1 

.5 

.VI/’ 

) 


' 





128 79 

I 

1 0 

10431 

1 

73.90 

18.56 

3*60 

21 15. S7 

ft » 0 

1 ’28 79 

1 

1 0 

10402 

1 

106*67 

38*35 

2190*53 

1954.14 

0.0 

128 79 

2 

1 0 

10402 ' 

2 

80.45 

22.30 

533.42 

207S.9S 

'3 • 0 

123 79 

2 

t 0 

13402 

2 

1 12.24 

43.32 

2036.59 

1967.21 

3.0 

123 79 

3 

1 0 

10402 

3 

94.32 

30.09 

844.90 

2052.45 

0*0 


?GET( INDEX# ' VGH- INDEX’ ) 

DD=32 NUMO* 3 CARD = 5 

C ‘3.305* 3.01/ 3.144. 0.35] 


? SUM (VGH. INDEX. SI » | ) 


*** KET 

SUM 

MEAN 

L0U 

HIGH 

STO 

3 

0 . 66039 3 E 06 

90.99 

->191*89 

126.53 

12.53 

9 

0.205569E 06 

25 * 72 

-3.24 

54 . id 

7.47 

10 

0.133710E OS 

1840.97 

-214. 89 

8517*24 ' 

1405*49 

1 1 

0.143352E 03 

1930.61 

1 ! .39 

2132*77 

100*96 

12 3.33 137 IE 03 

00=32 NUMO= 6 CAROs 
£ 8.553, 12.48/ 3 

0.05 

5 

.327. 43 

-1.53 

• 44 ] 

1 *89 

*3*35 


Figure 5—1 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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VGH DOMAINS 


1. 

1 

- 

Aircraft Type Identifier 

2. 

S 

- 

Aircraft Serial Number 

3. 

K 


A counter used to verify that data was 
input correctly into the information system. 

4. 

F 

- 

Flight Number 

5. 

G 

- 

Gust Indicator; 0 is no data, 1 is calm, 2 is rough. 

6. 

Cl 

- 

Clock - Upper 5 digits 

7. 

C2 

- 

Clock - Lower 5 digits 

8. 

IAS 


Indicated Airspeed (knots) 

9. 

Q 

- 

Impact Pressure (PSF) 

10. 

ALT 

- 

Pressure Altitude (Feet ) 

11. 

P 

- 

Static Pressure (PSF) 

12, 

ACC 

- 

Acceleration (g's) 


Figure 5-2 



In Figure 5-1 » we have a summary of the entire data hase, VGH, and 
we notice that we have some had data points. Under KEY=8 (IAS), the LOW 
figure is -191.89 and under KEY=10, the LOW figure is -214.89. Both of 
these fields should contain only non-negative numbers. Although it 
is possible to make corrections, it is somewhat cumbersome in STDS. As 
a result, we take the general attitude that bad data points will be 
deleted. In Figure 5-3, we isolate the negative IAS with the first RSLTA 
command. Since there is only one, we list the set for the purpose of 
examining it. Next, we isolate the negative altitudes. Since there are 
thirty, we only look at the first five. Next, with the relative complement 
command, RL, we form a new set VGH2 omitting the bad IAS data point. And 
again using the RL command, we form VGH3 omitting all the bad altitude data 
points. The PUT command allows us to save our updated data base as an MTS 
permanent file. Finally, we perform a statistical summary on the new data 
base in order to verify our deletions. 

Figure 5-4 describes housekeeping operations that may be performed 
from time to time. The LISTU command lists the sets that are currently 
active. Under the current version of STDS, the user may have only twenty 
sets active concurrently. As a result, it is necessary now and then to 
free sets that we are no longer using or save them as permanent files with 
the PUT command. The final GET command retrieves our updated version 
of the data base with our preferred name, VGH. 

Figure 5-5 describes the procedure for isolating items by serial 
number. Currently, the VGH Summary Program is used to obtain statistical 
information from the VGH data tapes. Each tape contains some flight 
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76C9MMENT ERROR DETECTION AND CORRECTION 


?RSLTAC8.VGH*0.0» I ASNEG) 

00=32 N'JMO= |2 CARD= 1 

C 0 • 320* 2*21/ 1 • 726* 6»88] 


7LISTC I ASNEG* 1 » l» • S \/' ) 

123 79 952 4 1 970253196 -191.39 


7RSLTAC 13»VGH»0.0, ALTNEG) 

00=32 NL'MDs 12 CARD- 30 

i 0*320* 0*32/ 1.707* 4.35] 


?LIST< ALTNEG* 1*5* * 


128 

79 

43 

3 

0 

128 

79 

l 19 

5 

0 

128 

79 

123 

5 

0 

123 

79 

126 

5 

0 

128 

79 

129 

5 

0 


l / ' > 


10402 

28 

79.72 

1 0002 

46 

76.14 

1 0002 

50 

72.28 

10002 

53 

62.43 

10002 

56 

68.73 


?RL( VGH* I ASNEG* VGH2 > 

00=32 NUMO= 12 CARO= 7262 
C 1.494* 6.76/ 1.869* 9.961 


3*24 3625.57 1853.51 


21.63 -186.11 2133.55 
19.71 -164.51 2128.39 
17.76 -214.39 2132.77 
13.27 -157.31 2123.34 
16.05 -99.63 2123.90 


7RLCVGH2* ALTNEG.V5H3) 

00=32 NUMDs 12 CARD = 7232 

E 0.360# 2.69/ 0.630* 5.66] 


-0.44 


0.0 

0.0 

0.0 

0.0 

0.0 


?PUT< VGH3* * NEWVGR * > 

% 00=32 NUMO= 12 CARD= 7232 
C 1 .590* 9.70/ 1 .965* 


11.42] 


7SUMCVGH3* INOEX* S2* 1 ) 
*** KET SUM 

3 , 0.6591 2 3E 06 

9 0.203090E 06 

13 0.133703E 08 

1 1 0.143196E 03 

12 4. 37 7 032 E 03 

oo=3 2 NU'iOs 6 Card= 
C 7.323* g.63 ' 


MEAN 

LON 

91.14 

30.43 

23.7 7 

3.14 

1843.34 

1 .28 

1 930.03 

1 1 .39 

0.05 

-1-53 

5 


524* 44.52] 



HIG3 
126.53 
54.73 
3517.24 
2116.12 
1 .89 


STO ♦ti- 
ll .92 
7.41 
1402.6 7 
1 00.74 
0.35 


Figure 5-3 


6 0 



?>C3MMENT HOUSEKEEPING OPERATIONS 
7LISTU 


NAME 

TTPE 

CRE* 

CAR)) 

01) 

NOD 

STATUS 

SIZE 

LOG 

VGH 

TAU 

1 

7263 

32 

12 

1 

35.11 

HdPTiVGH-SET 

INDEX 

'TAU 

2 

5 

32 

3 

r 

3.01 

VGH- INDEX 

SI 

TAU 

3 

5 

32 

6 

2 

0.03 

CORE 

I ASNEG 

TAU 

4 

I 

32 

12 

2 

3.01 

CORE 

ALTNEG 

TAU 

5 

30 

32 

12 

2 

0.35 

CORE 

VGH2 

TAU 

6 

7262 

32 

12 

2 

85.10 

CORE 

VGH3 

TAU 

7 

7232 

32 

12 

2 

84.75 

C0RE 

S3 

TAU 

3 

5 

32 

6 

2 

0.03 

C0RE 


?.FR£E(VGH> 
? FREE ( S I ) 


?FREF< IASNEG > 

?FREEE< ALTNEG > 

“FREEE ” IS -NOT DEFINED. 

7FREEC ALTNEG > 

?FREE( VGH2 > 

? FREE C VGH3 ) 

? FREEf S2 ) 

? GET i v' G T » * NE'.vV'Ga * > 

00=32 nu>io= i a card= 7232 

f 0 « 305» 0.03/ 0.362» 0.153 


Figure 5-4- 
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’^COMMENT ISOLATE 3* SERIAL NUMBER 

7RSEQC2.VGH# 13 5# SER 135) 

DD = 32 NUMi)= 12 CARO = 2500 

C 0.213# 1.37/ 1.719# 7.571 


? L I ST ( SER 135# 1#5# * /]/• ) 


ISO 185 

1 

7 2 

892516353 

78.64 

21-03 

1016-69 2339.63 

0.85 

150 185 

2 

7 2 

892516354 

79.01 

2 1 • 22 

1316-69 2339.63 

0.79 

153 185 

3 

7 2 

392516355 

79.9 1 

21.70 

1032.76 2038.45 

1 . 12 

150 185 

4 

7 2 

392516356 

63.31 

15.86 

1024.72 2339.04 

0.59 

153 185 

5 

7 2 

392516357 

70.44 

16.35 

961.34 2043.74 

0*76 

7 SUNK SER 135# INDEX# ! 

51 # 1 ) 



1 



*** XE/ 


SUM 


MEAN 


LOW 

HIGH 

STD ♦ 

8 

3.212446E 

06 

84.93 


60.60 

114.67 

6.16 

9 

0.617426E 

05 

24.70 


12.47 

44*9 4 

3*34 

10 

0.29251 1 E 

07 

1 170.04 


552.38 

2299.10 J 

319.59 

1 1 

3. 5070 IDE 

07 

2023.04 

1946.43 

2074.34 

30.17 

12 

0.273573E 

03 

3.1 1 


-1.53 

1 .89 

0.51 

00=32 NUMD= 

6 CARI) = 

5 






C 2.550# 2.60/ 2.751# 38.42} 


7RSEQC2# VGH#79# SER79) 

1)0 = 32 NUMD = 12 CARD = 2623 

C 0.227# 1.66/ 1.739# 6.53} 


LI ST ( SER79 

# 1 # 

5# 

’/!/•) 







128 

79 

1 

1 

0 

1 040 1 

1 

7 3.90 

18.56 

3-63 

2115.57 

0 * 0 

123 

79 

1 

1 

0 

l 0402 

1 

136*67 

38.85 

2190.53 

1954. 14 

0*0 

123 

79 

2 

l 

0 

10432 

2 

33.45 

22.00 

530.42 

2075.98 

0*0 

128 

79 

2 

1 

0 

13402 

2 

112.24 

43.02 

2306.59 

1967.21 

0*0 

126 

79 

3 

1 

0 

10432 

3 

94-02 

33.39 

844.90 

2052.45 

0*0 


?S'JM( SER79* 1NDEX#S2# 1 > 


!■> + !vE/ 

r* » 

JU * 


sAl? A M 

} IU' *1 ■« 

k 

U" r« 

H I GH 

ST9 

6 

0 • 25257 4E 

06 

96.11 

34.05 

126.53 

13.63 

9 

0.344214E 

35 

32.12 

3.94 

54.73 

3.44 

13 

0.493533E 

37 

1866.53 

l *28 

3517.24 

1559.04 

i 1 

0 « 524437 E 

37 

1973.33 

1 l .39 

2116.12 

116.52 

. 12 

3. 1 399 33E 

03 

3.04 

*■ v) • & <£ 

1 .02 

0.21 


00=32 NUMO= 6 CAR.) = 5 

C 2.665# 2.72/ 2.865# 33.60) 


Fig\ore 5—5 
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data for a particular serial number. The VGH Summary Program passes 
the tape and generates the summary. As a result all summaries are 
by serial number and it is difficult to obtain answers to queries that 
would require access to data for many serial numbers at one time. 

Appendix A gives some information that is included in the output 
record for the VGH Summary Program and we propose to answer a number 
of those types of queries using STDS; but in fact there is a large 
class of queries that would be useful that are not included in Appendix 
A because they cannot easily be obtained from records sorted by serial 
number. An example of one of these is: Give me all occurrences for all 

aircraft where the acceleration is greater than or equal to 1 g. This 
query is considered later in this section, 

However, in many cases it will be desirable to isolate data by 
serial number. First, we isolate the data for serial number 185 with the 
command: RSEQ(2,VGH,l85,SERl85) . Then, we list the first five elements 

and obtain a statistical summary. From this, we can directly obtain the 
average IAS (84.98) and the average altitude ( 1170.04) as well as the high 
and low values. We perform the same operations for serial number 79; an 
RSEQ followed by a LIST followed by a SUM. 

Figures 5-6 and 5-7 describe the procedure for isolating items by 
intervals in some domain. A number of items in Appendix A depend upon 
retrieval by interval; time interval, altitude interval and velocity 
interval. In this section, we demonstrate isolating readings by altitude 
intervals. The same technique may be applied for time and velocity 
intervals. 
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? ^COMMENT 


INTERVAL ISOLATION - ALTITUDE 


? 3RSAC l 0 « VGH* 2303 .0**3000 .0*1 NTERVAL ) 
00 = 32 NUMu = 12 CARD = 1243 

C 0*230* 1 *82 / 1.679* 6-855 


?L I STC INTERVAL* 

1 

» 5* 

* /{/' 

123 

79 

1 

1 

3 

13402 

123 

79 

2 

1 

0 

10432 

123 

79 

3 

1 

0 

13432 

128 

7 9 

4 

1 

0 

1 3432 

128 

79 

5 

1 

0 

13402 


t 136-67 

2 112.24 

3 98.34 

4 106.51 

5 103.17 


38.85 2193.53 
43.02 2006.59 
32.95 2322.57 
33.73 2113*52 
36*33 2134.51 


1954.J 4 

0.0 

1967.21 

0.0 

1966.33 

0 . 0 

1959.26 

0.0 

1958*13 

0.0 


7SUMC INTERVAL* IN0EX*S3> 1 ) 


*** ket Sum mean 

3 0 » I i 9562E 06 96.42 

9 3.401033E 05 32.34 

13 0. 3 02 50 4 E 07 2439.55 

1 1 0.239393E 07 1934.59 

12 0 . 323992E 31 0.00 

DD=32 NUMOr 6 CARD = 5 


C 1.266* 1.32/ 1.467* 37.203 


LOW 

34.81 

4.11 

2002.65 

11.39 

- 1.10 


HIGH 
121.58 
53 . 52 
2994.96 
1967.49 
J .44 


STD *** 
*13.80 
8.32 
293.33 
65.59 
0.23 


? INC INTERVAL* SER79* I NT79) 

DD=32 NlJMf) = 12 CARO* 514 
C 3*161* 0.23/ 0.294* 0.333 


7LISTC INT79 
123 79 

123 79 

128 79 

128 79 

128 79 


» 1 » 5* * /I /* ) 

1 1 3 13402 

2 I 0 10402 

3 1 9 13402 

4 1 3 10402 

5 1 0 10402 

♦ 


1 106*67 

2 112.24 

3 93.34 

.4 106.51 

5 133.17 


38.85 2190.53 
43.02 2936.59 
32.95 2022.57 
33.73 2118.52 
36-33 2134.51 


1954.14 

1967.21 

1966.03 

1959.26 

1958.12 


0.3 
0.3 
0.0 
0 . 9 
0.0 


7SUMC INT79* INDEX* S4» I ) 
*** KEY SUM 

8 0*5!3!90E 95 

9 0 * I 73330E 35 

10 •'*) * I 2 4 6 6 7 E 07 

It 9 • 993? 26 E 36 

12 M.592393E 01 

0D=32 NUM0= 6 C4R0 = 

C 3*527» 9.83/ 0 


MEAN L9tf 

O o , 3 a ^ /* . o \ 

34.70 4.11 

2425.44 2006.59 

1933.71 II.39 

0.01 -0.71 

5 

723* 35.291 


HIGH 

IOI Cl? 
• 4 - ov 

50.52 
2993.55 
I 967.21 
1 .00 


STD ■¥■¥* 

« /. /. “I 

* *1 • —I U 

8.9 7 
273.29 
3 7.12 
0.16 


Figure 5-6 
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?IN( £ NTERVAL, SER 1 8 5, I NT l 85 > 
00=32 NU’Ws 12 CARO* 7 1 


C 0.120 

9 

0. 

63/ 

0.239, 

0.771 





7LISTC I NT 135* 1 

.5 

,’/!/•) 






150 185 

1337 

3 

1 

970253353 

80.76 

22.13 

2126-21 

1 9 58 . 7 1 

0 • 40 

150 185 

1033 

3 

1 

970253351 

7 9,91 

21.70 

2159.15 

1956.37 

. -0.54 

150 185 

1089 

3 

1 

970253352 

85.41 

24.31 

2299. 13 

1946.43 

0.40 

1 53 185 

1090 

3 

1 

970253353 

85.73 

25.00 

2249.71 

1949.94 

- 0.41 

150 185 

1091 

3 

1 

970253354 

79.10 

21 .27 

2126.21 

1958-71 

- 0.56 


?SUmiNT185,INOEX,S5, l ) 


* KE7 

SUM 


MEAN 

L9tf 

HIGH 

STO 

8 

0 *6342735: 

04 

85.11 

63.25 

94.89 

5.45 

9 

0o 175677E 

04 

24.74 

13-60 

30.64 

3.04 

10 

0. 1 51323E 

06 

2133.36 

2032.65 

2299. 10 

79.88 

1 1 

0. 139006E 

06 

1957-33 

1946.43 

1967.49 

9.32 

12 

- .201997E 

01 

-0.03 

-0.76 

0.63 

0.47 


00=32 NUMD = 6 CARD= 5 

£ 3.079/. 0.03/ 0.279, 36.033 


Figure 5-7 
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First, we isolate all of the readings in the data base that have 
an altitude value between 2000.0 and 3000.0 feet. This is done 
with the BRSA command forming the set INTERVAL. We then list the 
first five items and perform a summary. In order to isolate the 
readings associated with serial number 79, we intersect the sets, 
INTERVAL and SER79 forming the new set INT79. This set contains all 
the readings for serial 79 in the altitude range 2000 to 3000 feet. 

From the cardinality indicator, we immediately know that there are 514 
readings in the set and a SUM command provides us with some additional 
information; the high, low and average IAS (121.58, 34,81, 99.84), and 
the high and low accelerations (—0.71, 1.00). One immediate advantage 
of an interactive system over a batch environment is the following. If 
the interrogator spots an area that needs further investigation, he has 
at his disposal the capability to immediately search to the root of 
the problem. 

Next, we isolate the readings for serial 185 by intersecting the 
set of interval readings, INTERVAL, with the set of readings for serial 
number 185, SERI85. As before, this isolates the readings for serial 185 
in the altitude interval . A LIST and a SUM provide a sample of the 
set contents and the usual statistical information. 

Figures 5-8 and 5-9 describe the procedure for isolating maneuver 
accelerations that exceed given thresholds. We will consider only the 
positive accelerations. However, we would normally be interested in the 
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? 5C0MMENT 


ISOLATE ACCELERATIONS 


?RSGTA( i 2 » VGM. 0 . 3. POSAC ) 

00=32 NU>1Ds 12 CARD= 1933 
C 3 »5 1 2. 3.92/ 2*000. 6*473 

? I N( P AS AC » SEP 1 3 5 » PCS AC 1 35 > 

DD=32 NUMO= 12 CARD = 1463 

C 0.246. 1.21/ 0.436. 1.563 


?SUMCpOSAClS5. IN0EX.S1. 1 ) 


*** KEY SUM mean 

8 0. 12284 l E 06 33-97 

9 0 -35303 1 E 35 24. 13 

10 0- 167533E 37 1145.17 

1 1 0.29693 IE 07 2029.94 

12 0 » 767 443E 03 0-52 

00=32 NUMO= 6 CARO= 5 


C 1.484. 1-55/ 1.685. 37.383 


LOW 

62.92 

13-45 

576-00 

1946-43 

0.39 


HIGH 
1 14.67 
44.94 
2299. 10 
2072.57 
1 *89 


STD *** 
6.69 
3*64 
309.32 
34.60 
0.15 


?IN(P0SAC.SER79.P0SAC79) 

00=32 NUMD= 12 C.4RD= 354 
C 0.151. 1-18/ 0.284. 1*373 


7SUMCPOSAC79. INDEX* S3. 1 ) 


*** k EY sum mean 

8 0.330032E 05 93.23 

9 0.106664E 05 30.13 

10 0.393555E 06 1 125.86 

35 0*71 9444£ 06 2032.33 

12 (3.17 1U48E 03 0.48 

00=32 NumO= 6 CARO= 5 


C 0.364. 0.4 4/ 0.566. 36.253 


LOW 
44.93 
6.84 
1 .23 
1836.25 
0.43 


. HIGH 
I 19.25 
43-57 
3875.88 
2116.12 
1 .02 


STO ■**■» 
12.33 
7.45 
316.17 
59.59 
0.1 1 


Figure 5-8 


6 7 



? c .COMMENT ISOLATE !G ACCELERATIONS (AND OVER) 

?RSGEA( 1 2»P0SAC» I .0*0NEG> 

DD=32 NUMD= 12 CARD= 34 

C 8.091* "0.10/ 0.193* g.g^i 


?L I ST ( ONEG » 1 * 5* ' /!/ 
123 79 1079 12 2 

128 79 II 54 21 2 

128 238 2089 40 2 

158 185 3 7 2 

150 135 7 7 2 


* ) 

852516353 112.17 
852516356 85*15 
852516356 72.44 
892516355 79.91 
892516359 74.47 


?S.UM(ONEG* INDEX* S3* | > . 


42.96 1402.94 2011,37. 
24.65 2864.72 1906.25 

1 7 .84 2117.13 1959.36 
21 .73 1032.76 2038.45 

18.84 937.76 2045.51 


1*02 
l .03 
l .44 
1 * 12 
I • 17 


*** KEY 

SUM 


MEAN 

3 

0*26061 5E 

04 

76.65 

9 

8 * 696438E 

03 

20.48 

10 

0.376145E 

05 

1106*3! 

1 1 

0*691 28 4 E 

05 

2033.19 

12. 

3*3331 93E 

82 

1*13 

0D=32 NUMDs 6 CARD- 

5 

l 0*04$ 

* 0*07/ 

3 

'•243* 36 


LOW 
64.58 
14.16 
764.32 
1906.25 
1 .00 


HIGH 
1 14.67 
44.94 
2864.72 
2053.45 
1 .89 


STD **+ 
12* 1 l 
‘7*34 
380.95 
27.73 
0.17 


7RSGTAC 12* ONEG * 1 «5*0NE5G) 
DD=32 NUMQ = 12 CARD= j 

C 0.004* 0.00/ 0.039, 

?LI ST ( 0NE5G » 1 *]*•/!/■ ) 

ISO 185 30 7 2 892516383 


0.093 

100.34 34.29 1024.72 2039.04 


1 *39 


Figure 5-9 
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magnitude of the a.ccelera/tion in both positive and negative directions. 

The positive accelerations may be isolated using the RSGTA command. To 
isolate accelerations for serial 185 we intersect the positive accelera- 
tion set, POSAC, with the set of readings for serial 185, SERI85. Statis- 
tics for this set may be obtained with the SUM command. Likewise, the 
accelerations for serial 79 may be obtained by intersecting the sets, POSAC 
and SER 79 and statistics may be obtained using the SUM command, 
may be obtained using the SUM command. 

As mentioned previously, it might be desirable to isolate all 
occurrences for all aircraft of accelerations equal to or over 1 g. This 
may be accomplished by using an RSGEA command forming a new set, ONEG. 

We know immediately from the cardinality indicator that there are 34 
occurrences. We may want to list all of them out at the terminal or 
only a few (as is shown) and perform a summary. The summary 
indicates that the high is 1.89; so we might decide to look at 
all over 1.5 g. This is accomplished using an RSGTA command and we 
find that the I.89 reading is the only one and we look at it using 
the LIST command. It occurred on serial 185. 

By now the reader will he familiar with the routine of (l) 
isolating a set describing a condition, ( 2 ) intersecting it with a 
set of readings for a particular serial number, and ( 3 ) obtaining 
statistics for the derived set. In this section, we follow an identical 
strategy in examining the rough air readings as shown in Figure 5 - 10 . 

Several items in Appendix A call for statistical analysis of 
sets that we have previously derived under Rough Air, Intervals, and 
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Accelerations. Figure 5-11 demonstrates the procedure for combining the rough 
air and altitude interval sets to determine the occurrences of rough air 
encounters in the altitude Interval 2000 to 3000 feet for serials 185 and 
79. 

The time and cost of this session is shown in Figure 5-12. Note that 
the connect time from initial dial-up to signoff was one hour and fifteen 
minutes and the central processor time was about 70 seconds on an IBM 
360/67. The dollar values are for a university environment and would 
be somewhat higher in a commercial environment. 
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? 'COMMENT 


ROUGH AIR 


?RSEQCS»VGH» 2 »R 0 UGH> 

00=32 NUMU= 12 CARD = 493 

t 0 . 164 * 0 * 58 / 1 * 596 * 4*933 • 


? !NCR 0 UGH»SER 79 *RUF 79 ) 

00=32 NUMO= 12 CARD= 146 
C 0 * 133 * 0 . 91 / 0 * 254 * 1*113 


?LI ST < RUF79 * 1 * 5* ' / 1 / ' > 


123 

79 

1 064 

3 

2 

892516353 

87*13 

25*86 

1967.61 

1970.35 

0*59 

128 

79 

1065 

5 

2 

852516353 

105*77 

33*21 

650 * 10 

2367.33 

3.53 

128 

79 

1066 

5 

2 

852516354 

96*95 

32.01 

807.41 

2055.26 

3.43 

1 28 

79 

1067 

5 

2 

852516355 

92*92 

29.40 

717.43 

2061 .99 

0.67 

123 

79 

1068 

5 

2 

852516356 

93*69 

33. 18 

754.94 

2059.19 

0.59 


?SUM(RUF 79 . INDEXaSl » 1 ) 


*** KET 

SUM 

~ 

MEAN 

LOW 

HIGH 

STO ?** 

8 

3 . 126637 E 

35 

86-72 

34.05 

1 19.25 

20.21 

9 

0 « 394350 E 

04 

27.01 

3.94 

48.57 

10.35 

10 

0 . 240504 E 

06 

1647.29 

8.60 

3875.38 

911.22 

1 1 

0 . 291 163 E 

06 

1994.27 

1836.25 

2115-57 

66.10 

12 

0 • 566793 E 

02 

0.39 

" 0.72 

1 .02 

0.40 


00=32 NUMDs 6 CARD= 5 

C 0 * 155 * 0 . 18 / 0 * 356 * 36*033 


? IN( RfiUGH* SER 185 * RUF 135 ) 
00=32 NUM 0 = 12 CARD= 301 


C 0 

-052 * 


0 . 

14 / 

0 . 375 * 

0 . 28 ) 


• 



LI ST C RUF 5 85 * 

1*5 


' - 





153 

185 

1 

7 

2 

892516353 

73*64 

21 .03 

1016-69 

2039*63 

0.8 5 

150 

185 

2 

7 

2 

892516354 

79*01 

2 l . 22 

1016*69 

2039*63 

0.79 

150 

185 

3 

7 

2 

892516355 

79.91 

21*70 

1032*76 

2033*45 

1 .12 

150 

185 

4 

7 

2 

892516356 

68.31 

15.86 

1024*72 

2039*04 

0.59 

150 

185 

5 

7 

2 

392516357 

70.44 

16.85 

961 *34 

2043*74 

0.76 


? SUMC R'JF 18 5 * INDEX *S 2 * 1 ) 


* * * KEf SUM 

8 0.226744E 05 

9 0 • 584959E 04 

10 0.237970E 35 

1 1 3.61 527 5E 36 

12 D. 21 4323 E 03 

00=32 NtJMD= 6 CAR 0 = 

C 0*313* 0*63/ 0 


MEAN 

75*33 

19.43 

956.71 

2044.10 

0.71 

5 


L 0 W 

62.92 

13.45 

533-87 

2017.90 

0.40 


516 * 36 . 55 ) 


HIGH 
I 14.67 
44.94 
1313.63 
2 v )7 1 .98 
1 .39 


STD * + + 
6.47 

3.52 

123.84 

3*53 

a . 2 1 


Figure 5-10 
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7SCQMMENT COMBINATIONS - ROUGH AIR IN THE INTERVAL 2-3 THOUSAND FT* 

?IN<RUF185» INT185* IR18S> 

NULL SET! 

C 0.022* 0-05' 0.113* 0.181 

? I N( RUF7 9* INT7 9. 1R79> 


DL>=32 

NUMD = 

12 CARO 

= 21 






C 0. 

033* 

0.20/ 

0.149* 

0.31] 



• 


7LISTCIR79* t » 

5* * / 1 / * 

) 






128 

79 1077 

9 2 

852516353 

44.15 

6*62 

2036-53 

1961* S3 

— 0.53 

128 

79 1078 

9 2 

852516354 

42.32 

6 . 00. 

2543*22 

1929.09 

-0*55 

128 

79 1118 

12 2 

352516392 

76.90 

20.11 

2222.55 

1951 .37 

0.79 

123 

79 1119 

12 2 

852516393 

40*47 

5.56 

2535.19 

1929.66 

-0*55 

128 

79 1120 

12 2 

852516394 

58.53 

1 1*63 

2495.06 

1932.51 

0.45 


? SUMC IR79* INDEX* S3* 1 ) 
**+ KE T SUM 

3 0.148883E 04 

9 0.395909E 03 

14 0 • 545536E 35 

11 0.404294E 05 

12 0. 560303 E 0» 
D0 = 32 NUMD= 6 C4R0= 

C 0.023* 3.33/ 0 


MEAN LOVJ 

70.90 34.31 

18.35 4.11 

2597*79 2036.53 

1925.21 1897.67 

0.27 -0.62 

5 

229* 36.03] 


HIGH STU +++ 

103.09 22.65 

36*27 10.90 

2985.50 ' 264.94 

1961 "53 13*35 

1 a 0 0 0 * 52 


Figure 5—11 
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?SC©MMENT END ©E STDS SESSION 


7S/STEM 

#UN 


9 S I G 

# ©FF AT 31:54.49 
9 ELAPSED TIME 
#CPU TIME USED 
#CPU STOP VMI 
#WAIT STAR VMI 
#L>RUM READS 
ft APPROX « COST Or 
#01 SR STORAGE- 
9 APPROX » REMAINING 


39-34-74 
74.383 MIN. 
69.983 SEC. 
239.933 PAGE-MIN. 
199*5^1 PAGE-HR. 


366 

THIS Run IS T.22 

3452.935 PAGE-HR • 
BALANCE: 5131 


Figure >-12 


S3 • 63 
S6 * 32 
513.00 


95 

51.14 

30 
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5.1 USING STDS FROM FORTRAN 


The reader will notice in Appendix A that there is a class of queries 
that are not directly answerable from an interactive system. For example, 
true airspeed does not appear explicitly in the data base; it must be 
calculated from impact pressure (Q) pressure altitude (ALT) and static 
pressure (P). The following equations in FORTRAN are used to obtain TAS. 
AARG = 518,688 - .OO35661T »ALT 
SS = 29.041+25 * SORT (AARG) 

XM = SQRT(5*(Q/PS + 1.0)**. 2857-1)) 

TAS = SS*XM 

In this subsection, we demonstrate how one may call STDS from an executing 
(FORTRAN) program. This provides the user with the capability to perform 
calculations based on items in his data base. Specifically, we demon- 
strate here a program that calculates average TAS in the altitude 
interval 8000. 0^9000. 0 feet. The interaction described herein was a 
session in itself from signon to signoff. The time and costs are given 
in Figure 5 - 15. 

The program was keyed in and debugged in previous sessions not detailed 
in this report. Figure 5-13 begins immediately after sign-on as we call 
the MTS text editor and request it to access the file TAS in which our 
FORTRAN code is stored. The command P *F *L is an editor command that 
requests that the text be printed from the first line (*F) to the last 
line (*L). The editor prints the entire program and finally is asked 
to terminate execution with the command, STOP. Now, we may consider the 
logic of the FORTRAN program itself. We assume the reader is familiar 



m TAS 

:P *F *L 

: X 

: 2 

: 2.25 

3 

: 3.25 

: 3.5 

: 3.6 

: 4 

: 4.1 

: 4.25 

: 4.5 

5 

: 6 

: 7 

: 7.25 

: 8 

: 9 

: 10 

: 11 

: 12 

: 13 

: 14 

: 14.25 

: 15 

: 16 

: 17 

: 18 

: 19 

: 20 

: 21 

: 22 

: 23 

: 24 

: 25 

: 26 

: 27 

: 28 

: 29 

:ST0P 


C** PROGRAM TO CALCULATE TAS 
INTEGER CARD 
CALL UN IV ( ’QQQQ ») 

• CALL GET(VGH, NEWVGH ’) 

IC=CARD(VGH) 

PRINT 94, IC - 

94 FORMAT (IX, ’CARD OF NEWVGH = * ,16) 

CALL BRSA { 10 , VGH ,8000 . 0 , 9000 . 0 , I NT ) 

IC = CARD (I NT) 

PRINT 95, IC 

95 FORMAT ( ’ CARD OF I NT = *, 16) 

TSUM = 0.0 

N = 0 

100 CONTINUE 
K = CARD(INT) 

IF( K .EQ. 0) GO TO 500 
N = N + 1 

C** GET Q PS AND ALT 

CALL GETELM(1,INT,33,4,Q) 

CALL GETELM( 1 , INT , 37 ,4 »ALT ) 

CALL GETELM(1,INT,41,4,PS) 

C** CALCULATE TAS 

AARG - 518.688-.00356617*ALT 
88 = 29 .04425*SQRT (AARG) 

XM = SQRT ( 5* ( (Q/PS + 1.0)** .2857 -1)) 

TAS = SS*XM 
TSUM = TSUM + TAS 

C** DELETE THE ELEMENT FROM THE SET 
CALL SUBSET(INT, 1,1,1, UNIT) 

CALL RL(INT,UNIT,INT) 

GO TO 100 
500 CONTINUE 
C** FINISHED WITH SET 
TMEAN = TSUM/N 
PRINT 90, TMEAN 

90 F0RMAT(1X, ’AVERAGE TAS = ’ ,F7.2,1X, ’KNOTS’) 
CALL EXIT 
END 


Figure 5-13 
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with FORTRAN and we will explain only the calls to STDS routines. At 
line 2.25, CALL UNIV ('QQQQ') sets up a universe named , QQQQ t . At 
line 3, CALL GET (VGH, 1 NEWVGH ' ) retrieves the data base 'NEWVGH' 
and assigns a local variable VGH for future reference. The cardinality 
of the set VGH is obtained with the CARD function and then printed by FORMAT 
9k. At line it, we make our primary call to STDS. We request that elements 
having altitude (domain 10) between 8000.0 and 9000.0 feet be placed in 
a new set called INT. The cardinality of INT is obtained and printed. 
Beginning at statement 5, we calculate average TAS for the elements of 
the set INT . FORTRAN statement 100 is the top of the loop that traverses 
the set INT, Statement 500 is the exit condition executed when the set 
INT is exhausted. 

The loop begins by obtaining the cardinality of INT. If it is non- 
zero we increment the element counter, N, Then we get the actual values 
for Q, ALT, and PS in the first element of the set by calling GETELM. 
Parameters three and four in the call are byte-oriented rather than 
domain-oriented. This is a human factors disadvantage of a number of 
the STDS functions. Then we calculate TAS and add the result to the 
total TAS, TSUM. Then we delete the first element of INT in the follow- 
ing manner. We isolate the first element of the set in the new set 
UNIT. Then we "subtract" UNIT from INT thus creating INT with one 
less element. At this point we continue by Jumping to the top of the 
loop. At statement 500, we have completed the task and we print the 
average TAS. 
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Figure 5-14 demonstrates how one compiles and executes the program 
from MTS, The command RUN -L+ RBLG: STDS SUB* + OLD: LIBRARY loads 
the file -L. In addition, it loads the STDS subroutines from the file 
STDSSUB* under the user ID "RBLG", Furthermore, those routines require 
some additional routines that are found on and loaded from the file 
LIBRARY under the ID "OLD". Details of this process may be found in MTS, 
Volume I [3 ] . 

After loading MTS issues the message "EXECUTION BEGINS". At this 
point our program is in control. We print the cardinality of the data 
base NEWVGH as requested in line 3.5 of the program. We print the cardinality 
of INT as requested in line 4.25. Finally, we print the average TAS 
for the altitude interval 8000,0-9000.0 feet as requested in line 2 6 
of the program. Figure 5-15 gives the time and costs for the session. 

The total cost of the session was less than $5.00 (at university rates) 
and the duration of the session was less than five minutes. 
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#RUN *FTN PAR=S=TAS ,L=-L 
# EXECUTION BEGINS 
NNO ERRORS IN MAIN 


#RUN -L+RBLG : STDSSUB*+OLD : L I BRARY 
#EXECUTION BEGINS 
CARD OF NEWVGH = 7232 

CARD OF INT « 39 

AVERAGE TAS = 118.23 KNOTS 

#EXECUTION TERMINATED 


Figure 5-14 
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#SlGNOl'F 


#Ohf A I (51:44.05 
#ELAPSLU TIME 
#CPU II ML USED 
#CPU STOR VMI 
#WAIT Si OR VMI 
#DRUM READS 


10-19-74 
4.754 MIN. 
29.592 SEC. 
32.867 PAGE-MIN. 
1.336 PAGE-HR. 


94 


# APPROX . COST OF THIS RUN IS 


#DI$K SIORAGE 


19.25 


# APPROX. REMAINING BALANCE: 


$4.68 

PAGE-HR. 

$257.28 


$.23 

$2.67 

$1.78 


$.01 


titiXBOIAL PAGE IS 
Of POOR QUALITY 


Figure 5-15 
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6.0 Concluding Remarks 


An information system incorporating the set-theoretic approach 
(e.g. STDS) provides a flexible environment within which to solve 
scientific problems needing data base support. The scenarios 
described herein suggest that the engineer or scientist can mani- 
pulate his own data base and obtain reliable answers to his queries 
in a simple and direct fashion. Furthermore, he may access his 
data base(s) from either executing programs or from on-line inter- 
active terminals and he may access his data base(s) using named 
variables rather than storage addresses. In addition, STDS provides 
a general capability to produce statistical information. 
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Appendix A 


VGH SUMMARY INFORMATION 


1. Aircraft serial No. 

2. Average tru airspeed over entire A/C flight history (total), 

3. Total flight miles flown, 

4. Total flight hours flown. 

5. Number of flights. 

6. Average true airspeed during check flight. 

7. Average true airspeed during operational flight. 

8. Flight miles during check flight. 

9. Flight miles during operational flight. 

10. Flight hours during check flight. 

11. Flight hours during operational flight. 

12. Number of check flights 

13. Average true airspeed during each flight condition (climb, cruise 
descent or operational - x-country) . 

14. Flight miles during each flight condition, 

15. Flight hours during each flight condition. 

16. Average TAS during altitude intervals. 

17. Flight miles during altitude intervals. 

18 . Flight hours during altitude intervals. 

19. Number of flights during altitude intervals. 

20. Percent of flights in time duration intervals (15 min. intervals) 

21. Average indicated airspeed - total. 

22. Average IAS during check flight. 

23. Average IAS during operational flight. 
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2b. Average IAS during each flight condition. 

25. Average IAS in rough and smooth air. 

26 . Average IAS in altitude intervals. 

27. Average IAS in type of air during altitude intervals. 

28. Percent of flight hours in velocity intervals (or 10, 20 know intervals). 

29. Total average altitude. 

30. Average Alt. during check fit. 

31. Average Alt. during operational fit. 

32. Average Alt. during each fit. condition. 

33. Average Alt. in rough and smooth air. 

34. Gust , pos. and neg# maneuvers, limit load factors. 

35. Occurrences of maneuver accels. in .1 g intervals. 

36 . Occurences of maneuver accels in check flight (.1 g intervals). 

37. Occurrences of maneuver accels. in operational flight (.1 g intervals). 

38 . Occurrences of maneuver Accels. in each flight cond. (.1 g intervals). 

39 . Occurrences of maneuver accels. in each alt. interval. 

40. Occurrences of gust accels. in 0 1 g intervals. 

41. Occurrences of gust accels. in each fit. condition. 

42. Occurrences of gust accels. in each alt. interval. 

43* Occurrences of derived gust velocities (4 FPS intervals). 

44. Occurrences of derived gust velocities in each fit. condition. 

45* Occurrences of derived gust velocities in each alt. interval. 

46. Occurrences of gust accels. during check fit. 

47* Occurrences of derived gust velocities during check fit. 

48 . Reading threshold. 

49 * Occurrences of maneuver accels in each velocity interval (10 or 20 knots)* 
50. Occurrences of gust accels in each velocity interval. 
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51. Occurrences of derived gust velocities in each velocity interval. 

52. Flight miles in each velocity interval, 

53. Flight hours in each velocity interval, 

54. True air speed in velocity interval. 
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